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ABSTRACT 



We give an exhaustive presentation of the semi-analytical approach to the model independent leptonic 
5^ I QED corrections to deep inelastic neutral current lepton-nucleon scattering. These corrections include 

photonic bremsstrahlung from and vertex corrections to the lepton current of the order 0{a) with 
soft photon exponentiation. A common treatment of these radiative corrections in several variables - 
leptonic, hadronic, mixed, Jaquet-Blondel variables - has been developed and double differential cross 
sections are calculated. In all sets of variables we use some structure functions, which depend on the 
hadronic variables and which do not have to be defined in the quark parton model. The remaining 
numerical integrations are twofold (for leptonic variables) or onefold (for all other variables). For 
the case of hadronic variables, all phase space integrals have been performed analytically. Numerical 
results are presented for a large kinematical range, covering fixed target as well as collider experiments 
at HERA or LEP(8)LHC, with a special emphasis on HERA physics. 
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1 Introduction 

The first deep inelastic scattering experiments have been performed at SLAC in 1968 and lead 
to the discovery of the parton structure of nucleons [1] . The cross section of the reaction of 
figure 1, 

eih)+pip,)^eik2)+Xip2), (1.1) 

was determined at a transferred momentum of \/Q^ ~ 3 GeV. The measurements were based 
on the registration of the scattered electron's energy and angle. The physical interpretation 
of the data took into account contributions from two important competing processes: the 
bremsstrahlung contribution to (1.1) of figure 2 with non-observed photon(s), 

e{h) + p{p^) -^ eih) + X{p2) + nj{k), (1.2) 

and also the elastic radiative tail, i.e. photonic bremsstrahlung in the elastic channel: 

e{ki) + p{pi) -^ e{k2) + p{p2) + n'-f{k). (1.3) 

Semi-analytical expressions for the cross sections of these radiative corrections to (1.1) had 
been calculated in the classical work of Mo and Tsai [2]. They were of considerable importance 
for the interpretation of the data. In the long period between the SLAC experiment and the 
start-up of HERA, the experimental devices for the study of deep inelastic lepton-nucleon scat- 
tering had a relatively simple structure. Therefore, there was little need to radically improve 
the treatment of the radiative corrections as developed in [2]. On the other hand, rising ex- 
perimental accuracy and considerably higher beam energies had to be met with corresponding 
improvements in the treatment of radiative corrections. 

Originally motivated by the needs of the BCDMS experiment [3], a series of papers [4]- 
[11] was published by the Dubna group in the seventies and eighties. They all are devoted 
to the model independent semi-analytical calculation of QED corrections arising from the 
reactions (1.2) and (1.3). The QED corrections connected with the electrons are treated for 
processes with virtual exchange of a photon between the lepton and the hadronic system. 
Applications to high energy muon scattering were first published in [5]- [7]. In [8], a technique 
was developed for a covariant treatment of the corrections, thus preventing the introduction of 
an unphysical parameter a); this parameter was used in [2] to divide the contribution from (1.2) 
into two parts which describe soft and hard photons separately. A covariant form of the peaking 
approximation of [2] and of soft photon exponentiation have been introduced in [9]. With the 
momentum transfer of \fQ^ ~ 10 GeV in the BCDMS experiment, the Z-boson exchange 
had to be taken into account in the Born cross section. Therefore, the structure functions 
became modified, and a new one describing parity violation had to be introduced [12]. All 
these developments, together with the contributions from QED lepton pair production [10] 
and the 0{a'^) corrections of the process (1.3) [11] have been included in the semi-analytical 
program TERAD86 [13]. A detailed comparison of the Fortran program FERRAD, which is based 
on the approach of Mo and Tsai with TERAD86 may be found in [14] . 

So far, the theoretical description of deep inelastic scattering relied on a picture, where 
leptons with some couplings (electrical charge, weak neutral couplings) interact via a neutral 
current with a hadronic system, being described by structure functions. These contain the 
quark couplings and distributions. In this sense, a factorization of the matrix element into a 
leptonic and a hadronic part occured. 



'{ki,m 




(/C2,m) 



p{pi,M) 




Xip2, Mh 



Figure 1: Deep inelastic scattering of electrons off protons: (a) Born diagram, (h) leptonic 
vertex correction. 



In a next step, the complete treatment of electroweak Standard Model corrections to 
deep inelastic lepton-nucleon scattering was performed in [15], thereby covering both the 
complete 0{a) photonic corrections with soft photon exponentiation and the full set of weak 
one loop insertions. Such a complete treatment of the weak corrections unavoidably spoils 
the aforementioned factorization between the electronic and the hadronic part of the cross 
section formula, which is the spiritual basis of the introduction of structure functions. See [16] 
for details. Therefore, the model independent approach to the radiative corrections has to be 
given up in favor of the quark parton model. The calculations in [15] go beyond the model 
independent approach in a second respect: A complete treatment of the photonic corrections 
includes photonic radiation from the hadronic state, containing as a part the electron-quark 
interference. The numerical calculations in the quark parton model approach were realized in 
the Fortran program ASYMETR [15]. At that time, the sophisticated treatment of the tiny effects 
which are covered in ASYMETR but not in TERAD86 was not justified by experimental needs so 
that the program ASYMETR did not find a broader attention. The same happened with the early 
study published in [17], where a leading logarithmic calculation of leptonic QED corrections 
for ep scattering at HERA energies was performed. Recently, initiated by the HERA Physics 
Workshops, the Dubna-Zeuthen group updated the treatment of the weak one loop corrections 
and the renormalization scheme used in ASYMETR. Further, the QED part was recalculated, 
carefully checked, and considerably compactified [16], [18]; the related programs are DISEPNC 
and DISEPCC. These two programs use the weak library DIZET [19]. 

Recently, the experimental techniques have been improved considerably and much higher 
beam energies are reached. In deep inelastic fixed target experiments, transferred momenta 
of v^ ~ 10 GeV and y/Q^ ~ 17 GeV were obtained at CERN [20] and Fermilab [21]. At 
HERA [22, 23, 24], the mass scale of the weak gauge bosons is in reach, \/Q^ ~ 100 GeV. 
As did the SLAG experiment, the recent fixed target experiments rely on the observation of 
the scattered electrons (muons) (with the notable exclusion of the neutral current neutrino 
scattering experiments [25]). By comparison, the HERA detectors HI [26] and ZEUS [27] 
represent a new generation. They allow to detect the scattered electrons, the hadronic final 
state, and also photons with high precision. Thus, for the physical analysis of the ep collisions 
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Figure 2: The two leptonic bremsstrahlung diagrams. 



at HERA one may use not only the familiar electron variables, 

pi{ki - k2) 



Qt = ih - k 



2) , 



yi 



Piki 



where 



-(^ki+pif = S + m'^ + M^ 



Xl 



AEeEp, 



01^ 
yiS' 



but also the kinematical variables from the hadron measurement, 



Ql = (P2 - Pi) 



Vh 



Pl{P2-Pl] 



Xh 



Ql 



Piki VhS 

or some composition of both, the so-called mixed variables [28, 30]: 



^m ^l • 



Vh, 



Xri 



VhS' 



:i.4) 



:i.5) 



:i.6) 



;i-7) 



Here E^, m and Ep, M are the energies and masses of incident electron and proton (see figure 1). 
Another useful set of hadronic variables has been introduced by Jaquet and Blondel [31]: 



QiB 



1 -1/h' 



yjB = Vh, 



xm 



VhS' 



1. 



Different choices of variables make no difference for the determination of the cross section 
of reaction (1.1) in the Born approximation. Although, there are huge differences in the 
predictions for the radiative corrections, because the kinematics becomes quite different. This 
may be seen from the tetrahedron of momenta which is shown in figure 3. For vanishing 
photon momentum /c, the simple Born kinematics is recovered. The differences concern the 
calculation of the corrections, but also, and maybe more important, their numerical values. 
In addition, the advanced HERA detectors allow the application of dedicated experimental 
cuts in the event selection which represent further potential problems for the calculation of 
radiative corrections. A satisfactory treatment relies almost exclusively on the use of the 
Monte-Carlo technique. Nevertheless, from the data with sometimes complicated cuts, it is 




P2 = Qh 



Figure 3: Spatial configuration of the momenta in reaction (1.2) in the proton rest system. 



usually possible to reconstruct some sufficiently inclusive and smooth cross sections. These 
cross sections then may be the subject of further study by a semi-analytical approach as is 
advocated here. 

The HERA physics workshops in 1987 [23] and 1991 [24] lead to an enhanced activity 
both in the calculation of radiative corrections and in the comparison between the results 
of different authors [32]. Several authors obtained new results for the radiative corrections 
using quite different techniques, which we quote here for completeness: leading logarithmic 
calculations [33]-[37], Monte-Carlo approaches [38]-[40], and semi-analytical approaches [16, 
18, 41, 42, 43, 44] for both the neutral current and the charged current reactions at HERA. 
Eventually, all groups were able to reproduce the numerical results of the above mentioned 
series of papers of the Dubna-Zeuthen group with reasonable precision [45] and, of course, to 
go beyond in many other aspects [32]. Considerable parts of the present article have been 
worked out during and after the 1991 workshop on HERA physics. 



In the following, we will restrict ourselves to neutral current scattering. The Born cross 
section is: 



|-.H^|:^,,,.,_i,e,,,,. 



;i.9) 



with the kinematical factors 






{y.Q' 



2(l-y)^^ 

2Q\2-y)S. 



The generalized structure functions Ai{x,Q'^) describe the electroweak interactions of the 
electrons with the protons via the exchange of a photon or Z boson and will be defined later. 
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It is well-known from the above mentioned earlier studies, e.g. from [16], that a treatment 
of the leptonic photonic corrections covers the bulk of the complete corrections to this cross 
section. Fortunately, both types of corrections - weak loop insertions and QED corrections 
related to the hadronic current - are relatively small. For a large part of the kinematical 
region they are below the experimental accuracy. If the experimental intention is a study of 
the hadronic current, one may be interested in a model independent description of the cross 
section which not necessarily uses the quark parton model. With the last remarks, the line of 
thought which we will follow from now on in this article is indicated: 

This article is devoted to complete, model independent, semi- analytical calculations of leptonic 
corrections to neutral current deep inelastic lepton-nucleon scattering in different kinematical 
variables. By complete we mean the full 0{a) corrections with soft photon exponentiation 
which are not restricted to the leading logarithmic approximation. By semi- analytical we 
understand that the Monte-Carlo technique is not used. We perform as many analytical 
integrations as is possible for a given set of kinematical variables. In order to get a double 
differential cross section, one has to perform three phase space integrations. In principle, 
one is interested to describe the cross sections with structure functions which may have an 
arbitrary dependence on the variables Xh, Q\- Then, for the case of leptonic variables, only 
one analytical integration, for the case of Jaquet-Blondel variables or for mixed variables - two 
integrations, and for hadronic variables all three phase space integrations may be performed 
analytically. 

In the above discussion, the characteristic elements of the calculation of real bremsstrahlung 
corrections have been introduced: 

• Choice of a reasonable phase space parameterization with a practical set of internal 
variables which are to be integrated over; 

• Choice of the order of integration and complete understanding of the corresponding 
kinematical boundaries, by necessity without neglect of masses; 

• Separation of the infrared singular part of the bremsstrahlung integral with use of a 
special rest system which has to be chosen appropriately; 

• Dedicated performance of the various hard bremsstrahlung integrations which are regu- 
lated or finite in the soft photon part of the phase space; 

• Calculation of the infrared divergent correction with a reasonable regularization proce- 
dure. Consecutively, compensation of the infrared singularity with that of the virtual 
corrections and elimination of the soft and hard photon separation with establishment 
of the lorentz invariance of the net correction. 



In order to give an impression of the spirit of the approach we quote here an expres- 
sion which contains the complete leptonic QED corrections to order 0{a) with soft-photon 
exponentiation to the Born cross section (1.9): 



d^a d'^aB 



dyhdQi dyhdQ 



exp 



-SinfiVh^Ql) 



2a 



3 3 



+ ^ E T^M^h, Ql)Siiyh, Ql) 
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The soft photon corrections factorize with the Born cross section. They are exponentiated in 
order to take into account the multiple soft photon emission: 









:i.io) 



The hard bremsstrahlung corrections are taken into account to order 0{a). They do not 
factorize but lead to modifications of the functions Sf while the generalized structure functions 
Aii^x^Q"^) remain unchanged: 



SiiVh^Ql) 



Ql 



^In^y^- 2^"^^^^ -Vh) -lni//iln(l-y/,) - -Li2(i//,) + -Li2(l) 



1 1 / 2 \ 
--ln?//,Lh + - 1 + — Lhi 

2 AX Vh 



Vhl ^' V ^Vh, 



S-2{yh,Q 



S^ --yh\Y?yh- (1 -i/h)ln^(l - yn) - 2{l - yh)\nyh\n{l - yh) - yhU2{l) 



(2 - ?>yh)U2{yh) + yhhiyhh^ + —{I - yh)Un - y(2 - yh) Iny/^ 



S^{yh,Ql 



SQl\-{2-yr.) 



--\v?yh + \n^{l-yh) + 2\nyh\n{l - yh) + 3Li2(t//, 



Li2(l) +hiyh + \nyhLY 



3 7 1 

-l/hLhi +yh- - + 2 (1 - 2yh) Int/^ L 



'1.111 



with Lh = \n{Ql/m^) and Lhj = U + ln[t//,/(l - yh)]- 

Since the structure functions depend on the hadronic variables and are hardly to be inte- 
grated over together with radiative corrections, such a compact result may be obtained only 
in hadronic variables. 

The purpose of the present article is fourfold: 

• We give a systematic presentation of the model independent approach to the leptonic 
QED corrections in terms of leptonic variables, thereby compactifying known results and 
transforming them into a modern terminology. 

• The complicated kinematics of deep inelastic scattering is explained in great detail for 
leptonic, mixed, hadronic, and Jaquet-Blondel variables. We systematically retain both 
the electron and the proton masses if needed. This is necessary for many intermediate 
steps of the calculation of QED corrections. 

• The complete leptonic QED corrections to order 0{a) with soft photon exponentiation 
in terms of mixed, hadronic, and Jaquet-Blondel variables are reviewed here for the first 
time^. 



^With the exception of references [46] and [47], which contain short collections of the main formulae. 
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• The contents of this article is the theoretical framework for the Fortran program TERAD91 
[41]. This program is used here for a systematic numerical study of the QED corrections 
for a wide kinematical range: fixed target deep inelastic /ip scattering, ep scattering at 
HERA, and ep scattering at LEP(S)LHC [48]. Naturally, the main emphasis will be on 
applications to HERA physics. 

The article is organized as follows. After an introduction of notations and of the Born 
cross section the radiative process (1.2) is described in section 2. In section 3, the phase 
space is parameterized such that the analytical integrations we aim at may be performed. 
A first integration will be carried out there. After this, the further treatment of the phase 
space will depend on the choice of kinematical variables. In section 4, we isolate and handle 
the infrared singularity, thereby defining soft and hard cross section parts in a covariant way. 
The soft photon exponentiation is performed. All this will be done separately for leptonic, 
mixed, and hadronic variables. The integrated cross section in the three above-mentioned 
variables is derived in section 5. There, for the case of mixed and hadronic variables, a 
second analytical integration is performed. Numerical results for the QED corrections at 
three different accelerator scenarios will be compared. In section 6, we introduce the necessary 
modifications for an accurate treatment of the photoproduction region. Section 7 contains a 
parameterization of the phase space, which deviates from that introduced in section 3. This 
parameterization will allow us to derive compact formulae for the QED corrections in terms 
of Jaquet-Blondel and hadronic variables. For the case of hadronic variables we perform all 
three phase space integrations analytically, thus obtaining a double differential cross section 
formula for the QED corrections without any integration being left. Section 8 contains a 
discussion. It includes a comparison with estimates which are based on the leading logarithmic 
approximation, a study of the influence of different choices of structure functions on the 
relative size of the radiative corrections, an example of cuts on the photon kinematics, and 
an outlook. In appendix A, we derive several phase space parameterizations, which are used 
in the calculations. In appendix B, the kinematical boundaries for the different variable sets 
are studied. The photons are treated totally inclusively in the main body of this article. 
Appendix B.1.2 is an exclusion to this. It contains some formulae with cuts on the photonic 
variables. Several tables of integrals which have to be used in order to perform the many 
integrations are listed in appendices D and E for the different phase space parameterizations. 
For completeness, we collect the cross sections in leading logarithmic approximation which 
have been used for comparisons in appendix F. 

2 Matrix elements and differential cross sections 

2.1 The Born process ep -^ eX 

In order to deflne the notations, we collect some basic formulae for the Born process (1.1). 
The matrix element which corresponds to the diagram shown in flgure la is 

s e^ 1 

M"" = J^{27r)''{p2\J^\p^)—u{k2){\Qe\l^ + XlA^e + ael5)}u{h), (2.1) 

where Ve and a^ are the vector and axial-vector couplings of the electron to the Z boson: 

Ve = l-4:\Q,\sm^ew, ae = l. (2.2) 
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Here, 6w is the weak mixing angle sin^ 6*^- = 1 — M^/M^ and Se|<5e| = Qe = —^ is the 
electron charge. Further, 

with the finestructure constant a = 1/137.06, the Fermi constant G^ = 1.16637- 10~^ GeV~^, 
and the Z boson mass Mz- In the numerical examples, we will use the following values: 
sin^^H- = 0.228, Mz = 91.173 GeV. 

The cross section of the reaction (1.1), 

das = 2'-^=^ Y. \^ ^n^^ (^i+Pi-^2-EpO, (2.4) 



V '"■-> spins 

depends on three variables, e.g. on S, the momentum transfer Q, and the energy transfer u: 
S = -2p,h, Q2 ^ (p^ _ p^y^ ^ ^ _2p^Q (2.6) 

In the following, we will use also the Bjorken scaling variables y and x, the parton momentum 
fraction, 

with 

Q^ = xyS. (2.8) 

The phase space element of the Born process is: 

_ dk2dp2,^^2,4 



dMt5\ki+pi-k2-p2)dVi 



2k^ 2p\ 

dydQ^dTh. (2.9) 



2^A 



In the last step, the integral over p2 has been performed with the aid of the (5-function. The 
integral over k2 is performed in appendix A.4.1. Further, we introduced the phase space 
element of the hadron system in the final state dT^ and its invariant mass M/^: 

^r. = n ^^'iP2 -Y.pd, Ml = -pI (2.10) 

In the following, the dV^ will become part of the definition of the hadronic structure functions. 
With our definitions we follow [49], 

§ E l^fl' = (^ {St:wl + 2x5«'Hi + x'5-H-;J) , (2,11) 
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where S^^ and VF^, A = 'j,I, Z, are the corresponding components of the leptonic and 
hadronic tensors. 

The phenomenological hadronic tensors W^^ [12], 






(2.12) 

describe the deep inelastic interactions of unpolarized nucleons with a photon and a Z boson. 
The real scalar functions W^ depend on the invariants Q^ and u. 

After neglecting in the contractions in (2.11) those terms from the squared Z exchange, 
which are proportional to alm'^, one observes a factorization of the leptonic and hadronic parts 
of the Born cross section: 

dVe ^^^'sJ:Mx,Q^)^^Sf{y,Q'). (2.13) 



The factorized functions Sf are: 

Sf{y,Q') = Q'-2m\ (2.14) 

Sf{y,Q') = 2[{1 -y)S'-M'Q% (2.15) 

Sf{y,Q^) = 2Q\2-y)S. (2.16) 

The factorized hadronic functions Ai{x, Q"^) describe the electroweak interactions of leptons 
with beam charge Qe and a lepton beam polarization ^ via the exchange of a photon or Z 
boson with unpolarized nucleons: 



A(a:,Q=^) = (2MWi) = 2^fC(^^g2)^ 
A(:r,g2) ^ j^ivW^) = ^^|'C(x,Q2), 



(2.17) 



The generalized structure functions J-'f'^[x,Q'^) are: 

^5^(a;,g') = Fi,2(x,g') + 2|Qe|(t;e + Aae)x(Q')G'i,2(a;,Q') 

+ 4 (vl + al + 2\v,ae) x\Q^)Hi,2{x, Q^), (2.18) 

T^i^^Q') = -2signiQi){\Q,\ia, + \v,)xiQ')xG,ix,Q^) 



+2 



2v,a, + X [vj + al)] x\Q^)xHs{x, Q^) I, (2.19) 



with A = C,sign{Qi). In the above formulae, the Qi is the charge of the lepton beam. 

The contribution of the weak loop corrections to (1.1) is not within the scope of the present 
article. We only mention that these corrections may be covered by an inclusion of real- valued 
weak form factors into the definitions of the weak neutral couplings (2.2). Such a program 
has been carried through first in [16], and was described and updated also in [50]. 
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The generalized structure functions J-'f^'^{x,Q'^) may also be used to describe some new 
phenomena, e.g. the virtual exchange of an additional heavy gauge boson Z' [51]. 

The running of the QED coupling a may be taken into account by a real form factor: 



a{Q' 



^Ff{Q' 



a 



l-EfQ}NfAFf{Q^y 



a 

IT 



5 Am 



9 + 3g^ + /n^ 



2my 



In 






Pf 



1 + 



4mj 



\^^ 



(2.20) 
(2.21) 

(2.22) 



The sum in (2.20) extends over all charged fermions /, Qf is the corresponding electric charge, 
and Nf the color factor: Nf = 3, 1 for quarks and leptons, respectively. The heavy fermions 
practically decouple in (2.20). For light fermions, with m^ << Q^ , the following approximate 
formula is valid: 



AF;(g^ 



a 
3^ 



In 



Q' 



m, 



5' 
3 



(2.23) 



The numerical factors in (2.18)-(2.19) are chosen such that some often used definitions of 
the structure functions in the quark parton model are matched. Assuming here the Callan- 
Gross relation, it is 



2xGi{x,Q^) 
2xHi{x,Q^) 



F2{x,Q') 

G2ix,Q') 

H2{x,Q^) 

xG3{x,Q^) 

xHs{x,Q^) 



xEW[q{x,Q') + q{x,Q^)], 
q 

X E \Qq\ Vg [qix, Q2) + q(x, Q% 
q 

xElU + al)[q{x,Q') + q{x,Q% 



X E \Qq\ a-q [q{x, Q^) - q{x, Q^ 
q 

xT.\ Vqttq [q{x, Q^) - q{x, Q^] 



(2.24) 



These definitions should help the reader to find a link to other articles on the subject, which 
often prefer the use of a slightly different notation (remind that we use the definition aj = 1 
for all fermions). A comprehensive presentation of the basic formulae for the description of 
deep inelastic scattering may be found in [52]. 

The allowed region of the variables y and Q^ is derived in appendix B.2.1 and shown in 
figure 4. In the following, sometimes the exact expressions in the electron and the proton 
masses are needed. So, in the figures with kinematics we will not apply the ultra-relativistic 
approximation. 

In figure 5, we show the Born cross section (2.13), rewritten as a function of x and Q^: 



(PaB 
dxdQ"^ 



27ra^ 



2\l-y-xy- 



M2 



m 



xQ^ [ \ ^ ^ ^ 

+ y{2-y)J^^^{x,Q') 



J^^^ix, Q'') + 2xy' II - 2^1 ^f ^(x, Q 



(2.25) 
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Q^ 



Figure 4: Physical region of the variables y and Q^ . 



At small Q^ the neglect of the Z exchange is a good approximation: 



dPal 



2-KC? 



dxdQ^ xQ^ 



M 



m 



2(1-1/)- 2x1/-^ + 1 - 2— 1 + 4a;' 



■ M'\ y' 



S 



Q' 



Q^ l + R 



J- 2 [x^ Q j, 



(2.26) 



where R is the ratio of the cross sections with virtual exchange of longitudinal and transverse 
photons, respectively: 



-<-^)^e^(--f)JISy-^ 



This cross section formula becomes in the ultra-relativistic approximation: 

Y^:n{x,Q')-y^n{^,Q') 



d'^aZ 27rQ;^ 



where the notations 



dxdQ^ xQ^ 



T2ix,Q' 



J^]{x, g") - 2xJ^^{x, Q') 
2xJ^{x,Q^)R{x,Q^) 



and 



Yi = l±{l-yy 



(2.27) 



f2.28) 



(2.29) 



(2.30) 
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Figure 5: Born cross section (Pa^/dxdQ"^ for deep inelastic neutral current scattering as 
function of In Q^ with x as a parameter. 

are introduced. The cross section (2.25) becomes in the uhra-relativistic approximation: 

d'^cTB ItioF' 
dxdQ^ ~ xQ* 

Here, we additionally assume the validity of the Callan-Gross relation [53]: 

J^^C{x, Q2) = 2xJ^^'^{x, Q^). (2.32) 

In the numerical examples, we will cover three typical kinematical regions: that of fixed target 
muon scattering {S = 1000 GeV^), of the ep collider HERA (^ = 4 ■ 30 ■ 820 GeV^), and 
of the ep collider project LEP(8>LHC (S* = 4 ■ 50 ■ 7000 GeV^). The solid curves in figure 5 
are calculated for LEP(8)LHC. They start at Q^ = 5 GeV^, the approximate lower limit of 
applicability of the structure functions which are derived from the available data. For a given 
value of X, the highest value of Q^ depends on S. This boundary derives from (2.8), taken 
at y = 1, and is seen in the figure. Further, the actual value of S has an infiuence on the 
cross sections (2.31) via the functions Y±. This leads to different predictions for the fixed 
target and HERA kinematics. The stars (for fixed targets) and crosses (for HERA) in the 
figure indicate the corresponding cross section values at the maximal allowed Q^ for given x. 
The first curve which shows a prediction for the fixed target case, e.g., is that for x = 10~^. 
Smaller values of x are outside the kinematical range for Q^ > 5 GeV^. The infiuence of S 
on the predictions is minor (in the variables and scales, which are chosen here). In the double 
logarithmic representation, the cross sections are nearly linear functions over the complete 
kinematical region. The cross section behaviour is determined by the characteristic l/{xQ'^) 
dependence. The l/Q"^ arises from the photon propagator (and, at very high energies, also 
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from the Z propagator). The event rates are highest for the smallest values of x and Q^ . 
With the lower bound on Q^ being fixed, the accessible kinematical range in x and Q^ rises 
considerably with the available beam energies, although the event rates at the highest values 
of X and Q^ are substantially suppressed. 

In sum, the ep collider with the highest accessible beam energies will open the field not 
only for a study of very deep inelastic scattering, but also, at high rates, the study of the 
region of very small x. 

The cross section predictions may depend crucially on the choice of structure functions. 
For the illustrational purposes, we decided to use the MRS D'_ parameterization in the DIS 
scheme [54] as it is compiled in [55]. The influence of the structure functions on the theoretical 
predictions, including the numerical values of the QED corrections, will be a subject of the 
discussion. 



2.2 The radiative process ej9 -^ eX^ 

The differential cross section for the scattering of electrons off protons (1.2), originating from 
the Feynman diagrams of figure 2, may be derived from the following matrix element: 



M^ 



"'^''''-Mk){27.f{P2\JM) ^ 



(27r)V2 

X U{k2)\ \Qe 



Ql 
1 /_. . X 1 



■7^ \2kia - k-iaj [2k2a + 7a^j 7/^ 



+ X 



1 ^ 1 ^ 

—l^l {Ve + OeTs) i'^kia " ha) ('2k2a + lak) 7/. (^e + ^£75) 



>u{ki), 

(2.33) 



where the electron propagators have the denominators zi or Z2'- 

zx = —2kik, Z2 = —2k2k. 

The corresponding cross section^. 



(2.34) 



^(27r)^r)9 v-^ I . ,R 2 ,^, , ,^-^ ,.dk2 dk -n- dpi 

V^S spins i Z,K2^K • ZPi 



may be rewritten as follows: 



dav 



ia" 



TT'V^Qt 



SdT, 



where 



s = Im {s;^w;:^ + 2xslwl, + x's^w^ 



(2.35) 



(2.36) 



(2.37) 



^The squaring of the Feynman diagrams has been performed with the program for algebraic manipulations 
SCHOONSCHIP [56]. 
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and 



dT 



dk2 d'p2 dk 
2k^ 2p9, 2A;0 



dM^S^{k2 +P2 + k-ki-pi). 



(2.38) 



The sum (2.37) may be treated exactly in the masses of electron m and proton M. As 
mentioned above, the Born cross section (2.13) exhibits a nice factorization property. We 
have shown explicitly, that the same is true for the radiative cross section after a neglection 
of the tiny terms proportional to x^fle^^^/-2i(2) and x^Ogm^/zf(2)- 



S =S{£,I,Zi,Z2) =J2'^ii^h,Ql)Si{£,I,Zi,Z2). 

i=l 



(2.39) 



The neglected terms are due to the Z exchange. They are vanishing at small Q^ since the 
Z propagator is suppressed there compared to the photon propagator. So, the corresponding 
approximation is sufficiently well even in the photoproduction region. In that region one has 
to retain, however, all particle masses in the cross section contributions from photon exchange. 
We should mention in addition that certain bremsstrahlung corrections, which are proportional 
to m'^/z^2) yield finite contributions to the cross section even at larger Q^; this may be seen 
from an inspection of the tables of integrals in appendix D. 

Besides on zi and Z2, the radiator functions Si{S,X, zi, Z2) depend on two additional two- 
dimensional sets of variables S and X: 



Si{S,I, Zi, Z2) 



l + J_Wt-4,„^) + (i 



1 

Z2 



1 



{Ql + QJ) - 2m^ 



S2\£.,T., zi, Z2) 



2M' + -^{qIHi - yH)S' + (1 - m)il -m + y,)S' 
Z1Z2 I 



(2.40) 



Z1Z2 

r2/'/n2 I o*„2m I o,»„2 



2M'{Ql + 2m')] + 2m' 2(1 - yi) + y^iyi - yn) 



S' 



1 

Zl 



y^S' + M\Ql + Q't 



I 

Z2 



y,(l-y,)S'-M\Ql + Q't 



2m' ^ 

zi 



2m' ^ 

zi 



2/n2 

h 



:i-l/z)(l-l// + t//.)^'-M2Q 



[\-yy,)S^-M''Ql 



(2.41) 



Sz\£-,1, Zl, Z2) 



S\2QI {Ql + 2m^){2-yi) 



1 
+ — 

Zl 



u. 



Z2 



2m'Q 



Z1Z2 
2Ql-yH{Ql + Q' 

\i2-2yi+_ 

Zl Z2 



2{1 -y,)Ql + yf,{Ql + Qj 



2,02 



\i2-2yi + yh) + \{2-yh) 



(2.42) 



The first two of the functions Si are the result of an exact calculation in both the masses 
m, and M for the case of photon exchange. For Z exchange, they contain certain harmless 
approximations in the electron mass, which were mentioned before (2.39). The same is with 
the ^3. 
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These radiators Si will be the central objects of the derivations in the following sections. In 
their arguments, the symbol S stands for two external variables, which will not be integrated 
over, while X denotes two internal variables. These latter two become, together with Z2 (or 
Zi), part of the integration measure. In the first part of the paper, the variable sets S and X 
will be selected out of the following set of variables: 






(2.43) 



The radiators Si are the main ingredient of the subsequent integration over the radiative 
phase space. The functions Si and S2 agree analytically with expressions obtained earlier [6]). 
Further, the Monte Carlo program HERACLES [38] is based on formulae, which are equivalent 
to our functions Si (i=l,2,3) [57]. 

The following sections will be devoted to a number of analytical integrations, aiming at 
semi-analytical, compact expressions for the QED corrections. 

2.3 The elastic radiative tail ep -^ epj 

A possible background process for deep inelastic scattering is the elastic scattering (1.3) from 
which the so-called elastic radiative tail originates. The corresponding formulae have been 
derived in [5, 42]. They may be obtained from the formulae for reaction (1.2) applying the 
following relations between the generalized structure functions Ai{xh,Qh) (i=l,2,3) and the 
generalized elastic form factors Ai{Ql): 

Ai{xh,Ql) = xhAi{QI)6{1-xj^), (2.44) 

A2,3{xh,Ql) = ^A2AQI)S{1-xh), (2.45) 

where Ai{Ql) are defined by formulae (20)-(22) of [42]. 

Further, from (3.13)-(3.16) of the subsequent section, one may get formulae (26)-(28) 
of [42], the expressions for Sf2 3 in the ultra- relativistic approximation. The exact in m and 
M expressions for 5^2 niay be found in [5], equation (37). 

3 A covariant treatment of the phase space 

3.1 Phase space parameterization 

The cross section of reaction (1.2) is characterized by six independent invariants which may 
be taken as follows: 

S, Vh Qf, Vh, Ql, Z2. (3.1) 

In addition, we use 

zi = Z2 + QI-QI (3.2) 
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We will now rewrite the phase space of process (1.2) such that it is expressed by these 
variables. The phase space is defined as follows: 

where the dT^ is already split away. After rewriting d'p2/2p^ = d'^p2 ^{pI, + ^h) ^"^^ taking the 
integral over d'^p2 using the 5 function, we arrive at 



rlh rlh 

dV = ^,§^S[{p, + Q,y + M-^]dMl (3.4) 



With a trivial calculation, which is explained in appendix A.4.1, one gets 



dk2 ttS 2 

Zhjrf Z\/ A 



dyidQl (3.5) 



Introducing the notation 

dk 



dT, = ^6[ip, + Q,y + M^]5[Ql - (p^ - Pi)')], (3.6) 



we arrive at: 

ttS 



dV = — ^ dyi dQf dMl dQl dVu- (3.7) 



With M| = —pi = M^ + yhS — Ql, the differential dM^dQl may be replaced by SdyhdQ\. In 
appendix A. 4. 2, we derive the following parameterization, which will be used for the subsequent 
integrations: 

^r, = -^, (3.8) 



where Rz is a quadratic polynomial in Z2- 

In sum, the phase space in terms of the invariants yi^Q^^yu^Q'h, together with one addi- 
tional variable Z2 (or zi), is [6]: 

:.2 



TtS* f 2 1 7/-i2 dZl{2) 



f dmdQUyhdQl^. (3.9) 



3.2 An analytical integration 

Taking into account (2.39) and (3.9), the differential cross section (2.36) for the process (1.2) 
can be written as follows: 

da^='^-^S{E,Z,z^,Z2)dEdZ-^^^. (3.10) 

Here and henceforth we use the following variables: 

X = Ii = {yh,Ql); 

(3.11) 



leptonic variables : 


S = 


-- Si -- 


= {yi,Ql), 


X = 


Ti -- 


= iyh,Ql); 


mixed variables : 


S = 


~ ^m 


= (yh^Qf), 


X = 


1m -- 


= {yuQD; 


hadronic variables : 


8 = 


-- Sn -- 


= iyh,Ql), 


X = 


Th -- 


= {yuQD- 
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For the case of mixed variables, we define in addition: 



Jyr. 



VhS' 



(3.12) 



The physical regions of Sa and Ta, a = l,m,h, are determined in appendix B. Here, we would 
only like to mention that the variable Xm niay not be restricted to the interval [0,1]; see the 
discussion in appendix B.3. 

For the calculation of the double differential cross section (1.2) one has to perform a 
threefold integration over the squared matrix element (2.36) with the following integration 
variables: Z2 (or zi), and the two invariants in X^, a = l,m, h, respectively. 

The first integration is that over Z2. Since the generalized structure functions AilxhjQl) 
are independent of -2i(2), it is sufficient to integrate the kinematical functions Si{£,I,Zi,Z2) 
of (2.40)-(2.42). The integration limits may be found in (B.9), and the necessary table of 
integrals in appendix D.l. 

The following expressions are the result of the integration: 



S{S,I) = ^f^S{S,I,z„Z2)^J:A{xf,,Ql)S,{S,I) 
where the Si{S,X) are 

B2 1 1 



(3.13) 



S2{Ea) 



VC21 
1 



Ql-Ql 



2m" 



cr 



lA. 



yi) 



(3.14) 



^ 



+ 



M\Ql + Ql)-yn{l-yi)S' 



Ql{l)[{l)-yH]S' 



{Ql-Qf)VC-2 
+ (1 - m) [(1 - m) + y,] S' - 2M\Ql + 2m') 



2m'S 



2o2 



- 2m' 



B. 



[(i)-y,][(i-i/0 + i//.] + (i)(i-i//) 

2M' 



^3/2 



2/n2 



[\)\{\)-yu\S'-M'Q 



'2 



M; 



S,{£,I) 



S 



2Ql{Ql^2m')\{\)^{\-y{)\ 



(3.15) 



Ql-Ql 



2il-m)Ql-y,iQl + Ql 
- 2m'SQl^^ [2(1) - y,]| + |(1) - -(1 - I/O}, 



(3.16) 



where [6] 
A2 = X 



qi 



22 



2 ^2\2 



B, = [2M'Q^{Qf-Ql) + {l-yi){yiQl-y,Q^)S' + S'{l)Q^{yi-y,)} 

C2 = [{l-yi)Ql-Qn{l)-yk]fs' + 4m'[{yi-y,){mQl-y,Qf)S'-M'{Ql-Q, 
^ C^{{l)^-{l-yi)}. 

The kinematical function Xg and the coefficients Ai, Bi, Ci, i=l,2, are derived in appendix A. 
Again, for the photon exchange and the 'jZ interference cross sections, the radiators are exact 
in the masses m, M. 

Formulae (3.13)-(3.16) describe the contribution of the radiative process to neutral current 
deep inelastic scattering without applying a cut on the photon kinematics. In (3.13), an 
integration over the azimuthal angle ip of the emitted photon [as introduced in (A. 21)] is 
already performed. In principle, one could proceed here in a different way and apply some 
simple photonic cuts. This is done in section 8.3 and appendix B.1.2. 

Further, we should indicate that the resulting expressions are singular at /c — > 0. There, 
Zi and Z2 vanish. It may be seen from (3.2) which relates zi and Z2 that Q| and Qf become 
equal in this limit. The same happens with yh and yi [see (2.43)]. This reflects the infrared 
singularity of the radiative cross section which deserves a special treatment before we can 
perform additional numerical or analytical integrations. 

4 Removal of the infrared divergence 

4.1 The infrared divergence 

In the last section, it has been shown that the double differential cross section of the pro- 
cess (1.2) can be written in the following form: 



dan _ 2a^S^ 
He ~ Ac 



fdlJ2A{x,,Ql)^S,{S,I). (4.1) 



As was mentioned at the end of the foregoing section, the integrand of (4.1) diverges if the 
kinematics is such that the photon momentum vanishes. Then the kinematical functions 
Si{S,I) become infrared singular as may be seen from formulae (3.14)-(3.16). Usually, in 
complete 0{a) QED calculations the infrared singularity is treated as follows. One introduces 
an (artificial) infrared cut-off parameter e, which divides the photon phase space into two 
parts: In the soft photon part, the photon momentum may become infinitesimal, while in the 
other, the hard photon part, it is enforced to remain finite. The soft photon contribution to the 
cross section is treated analytically; it depends on the cut-off which is introduced in a specially 
chosen lorentz frame. Thereby the lorentz invariance is spoiled intermediary in the calculation. 
The soft photon contribution contains the infrared singularity, which is compensated by the 
vertex correction of figure lb. The hard photon contribution is calculated by some dedicated 
integration method and depends also on the cut-off e. Finally, the compensation of the various 
occurrences of the infinitesimal parameter e is performed numerically, often with the need of an 
adjustment of it to the kinematical situation. Examples of such an approach are [2] and [43]. 
For the calculations, which are aiming at the leading logarithmic approximation, the soft 
photon problem may be solved completely by an analytical calculation of all the contributions. 
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This implies both an exphcit analytical compensation of the infrared cut-off parameter e and 
the semi-analytical integration of the hard photon bremsstrahlung. 

Here, the infrared problem will be treated with a covariant method [8] with similar features 
as that used in the leading logarithmic calculations. As mentioned above, a cut-off parameter 
e will be invented. Then, from the exact squared matrix element the terms, which contain 
the infrared singularity are extracted. These are of much simpler structure than the exact 
expression. The resulting simplified cross section part will be integrated analytically over the 
full phase space, thereby compensating the infrared singularity with the vertex correction and, 
more complicated, eliminating the cut-off parameter e. By construction, the infrared divergent 
bremsstrahlung cross section part contains hard and soft photon contributions. This is quite 
similar to the leading logarithmic approximation. 

The rest of the bremsstrahlung contribution is free of the soft photon problem. The 
integrand is quite complicated. Nevertheless, one may perform analytical integrations over all 
those internal variables, which are not an argument of the structure functions. 

One may represent (4.1) in the following form: 

dS \ d£ d£ ) d£ d£ d£ ' ^ ' ' 

where da^ /d£ is the infrared divergent part to be defined below and da^/d£ is finite at A; — * 0. 
Such a separation is not unique. The essentials of the approach will be the same for all the 
different ways of integrations over sets X of invariants, although certain details of the specific 
calculations will differ. 

Going back to (2.36) and (3.7), the fully differential cross section may be written as 

d(TTx = — ^ — —j;S{£,X,zi,Z2)d£dIdTk. (4.3) 

From (2.40)-(2.42), one may get the limit of S{£,X, Zi, Z2)/Qt at A; ^ 0: 



lmi^5(^, J,zi,Z2) = JZA{x,Q^)^^Sf{y,Q'):F'\Q\z„z,), (4.4) 



r-togr^ ' ' " ^^ tr'^'^'Q'' 



where 



TR.^9 , Q2 + 2m2 „/l 1\ , , 

^'""(Q', ^1, Z2) = ^ - + - ' (4-5) 

Z1Z2 \Zi Z2 ) 

and X, y and Q^ are variables as defined in the Born kinematics. The factorized universal 
function JF^^ is the well-known Low factor [58], 

T^^=(^ ^V (4 6) 

\lkxk 2k2k) ' ^ ^ 

and contains the infrared singularities of the photonic bremsstrahlung from the electron line. 
Now, we define the infrared part of the cross section as follows: 

^ = ^ I dIdT,j:Mx,Q')^^Sfiy,Q')J^'^{Q^z^,Z2) 

d'^a-B 2aS 
d£ 7r2 

- '^-sfin (4.7) 

at 11 
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JdIdTkJ^'''{Q\z^,Z2) 



The integration range in (4.7) will be split into two parts: 

TX J 



+ _ / dldTkJ'^''{Q\z^,Z2)e{h - e) 

TT J 

Cft(^,e) + Cd(^,e). 



(4.8) 



Here, the ko is chosen to be the energy of the emitted photon in the proton rest frame (A. 4), 

S 



k' 



2M 



iVi -Vh), 



(4.9) 



and e is an infinitesimal parameter, e > 0, e -^ 0. The first term will constitute the soft photon 
contribution to (5soft(^)- 



Cft(^,e) 



2S 



2 r dk 

IT 



dA 

2ko 



{pi + Qhf + Ml 6 Ql-{P2-Pif ^'^(g^^l,^2)^(e-A;o 



j^^:F'\Q\z,,z,)e{e-k,). 



(4.10) 



The second part contains an integration over nearly the complete phase space and implies also 
contributions from hard photons, but with a considerably simplified integrand compared to 
the complete bremsstrahlung integral. It will also be dependent on the infrared cut-off e, but 
is infrared finite for finite e: 



"hard (^5 ^, 

With the additional notation 
^^^(^,J) = 



5" f ,^dZi{2) ^lRf^2 



dl ^ 

71 J ^JR 



T'''iQ\z,,Z2)eiko-e). 



(4.11) 



-/%^-^^W,-l,^2) 
TV J V-Kj- 



g2 + 2m2 / 1 



1 



Qi ~ Qh yCi yC2 



— m 



Bi 



C'l 



3/2 



-D2 ' 



(4.12) 



the hard part of the infrared divergent correction becomes: 

Cd(^,e) = JdIJ''''{S,I)e{ko-e). 



(4.13) 



The integration of (4.5) over Z2 has been performed with the table of integrals D.l. The 
coefficients Bi,Ci, i=l,2 are introduced in appendix A. 

For later use, the infrared divergent part of the corrections is quoted also in a slightly more 
general notation: 



C(^) 



TT^ 



s 



/f^^^(Q^^l,^2) 



(4.14) 



where the Q^ on the right hand side is, by definition, one of the external variables and the 
T^^iQ^,zi,Z2) is defined in (4.5). 
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4.2 The soft and hard parts of ^|^ 

The soft photon contribution 6l^f^^{S,e) will be considered first. The integration region is 
limited to photons with an energy, which is smaller than e. After inserting JF^^ as defined 
in (4.5) into (4.10), one gets 



Cft(^,e) = 2{2nrJ 



d^k 



(27r)3A;0 



-k' 



g2 + 2m2 



m 



m 



[-2kik){-2k2k) {-2kiky 



-2k2ky 



(4.15) 



The expression (4.15) will be regularized now by changing into the (n — l)-dimensional 
space [59]. After an integration over (n — 3) azimuthal angles, two integrations remain to 
be performed: one over the polar angle d of the photon and the other over the photon energy 



"soft(^) 



2(2 



TT 



da- 



, ^, ^, , , , ^^ , , (ky-^dk'^ / (sin^)"-^^^ 

(2v^)"r (n/2 - 1) Jo /i"-4 Jo^ ' Jo ^ ' 



g2 + 2m^ 



m 



m 



(A;0)2(l-/5„cos^«)2 (A;?)2(l-/3iCos^i)2 {k^^)^{l - p^cos^^y 



(4.16) 



In order to simplify the angular integration, in the first term a Feynman parameter integration 
is inserted, which linearizes the angular dependence of the denominator: 



da 



-2kik){-2k2k) Jo [{-2kik)a + {-2k2k){l - a)^ 



(4.17) 



Further, the relations 



ka = kia + k2{l — a). 



-kik = A;°A;°(l-Acos^i), 



A 



kr 



i = 1,2, a. 



(4.18) 

(4.19) 
(4.20) 



are used where the (3i are (n — 1) dimensional velocities, and the 'di the corresponding spatial 
angles between the photon three- momentum and ki. We now go into the rest system of the 
proton, where it is pi = 0. In this system, the j3i are expressed in terms of the invariants (A. 4). 
We get: 



/i"-4 Jo 



0n„0\n-5 



dU\k- 



i^/f^) 



n—4 



n- A 

1 



n-A 



1 + (n - 4) In - + 



(4.21) 



Soft photon emission is isotropic. For any of the three terms under the integral, one has the 
freedom to choose a coordinate frame in the proton rest system with the (n — 1)** axis being 
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parallel to ki, i.e. 'di = i). After a trivial change of variable, cos-i? = ^, 

r-TT (sin^)"-3rf^ /-l (1 - ^2)n/2-2^^ 



{1 - l3i COS i!}y 



-1 (1 - m' 



l + -(n-4)ln(l-a + - 



-i(l-AO' 
From the above substitutions, one obtains the following expression for 5^^^^_{£,e): 



(4.22) 
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(4.24) 



(4.25) 



where the e is a parameter, which is often used in the MS subtraction scheme. The V^^ 
contains the pole term, which corresponds to the infrared divergence, as well as the Euler 
constant 'Je- The arbitrary parameter /i has the dimension of a mass. 

After the integration over C, and a with the aid of the integrals of appendix D.2 one gets: 
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2e 



(Q2 + 2m^)L^ - 1 






with iS$ defined in (D.22) and (4.31). Further, 
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(4.28) 
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The equation (4.26) with 5$ defined in (D.22) is exact. In the ultra-relativistic limit it is 



(4.29) 
(4.30) 
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and one gets the following short expression: 



Cft(^,e) 



2e 



2pm ^ 21n In 

/i 



M2g2 
m2M2 



(4.31) 






+ iln2^- hn\l-y)-U,{l), 



m^ 



(4.32) 
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where 



■^0 y 



(4.33) 



The infrared pole V^^ in (5s^t(^, e) has to cancel against a corresponding contribution from 
the QED vertex correction: 
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(4.34) 



This cancellation may be seen explicitly from the above expressions. 

The parameter e in (4.32) has to be compensated by corresponding terms from 5^-^^{£, e). 

In order to obtain the 5^^^{£, e), one has to perform the integration of JF^^(£^,X) over the 
physical regions of the variable sets Xa{a = l,m,h) as derived in appendix B, but with the 
exclusion of a small region around the infrared point JF; see also the figures in that appendix. 
These integrations are tedious and will be commented on in appendix D.3. 

In the ultra-relativistic approximation, one gets the following expression for ^^^j-d i^ leptonic 
variables: 
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In mixed variables, it is: 
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(4.36) 



Finally, in hadronic variables: 
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(4.37) 



At this intermediate state of the calculation, we would only like to remark that the depen- 
dencies of (^Jj^rd on the ln(2e/?72), hi{S/irP')^ and \n{S/M'^) are fictitious. They are compen- 
sated by corresponding dependencies of (5^^^ completely; see (4.32). At the other hand, the 
\n{Q'^ /irP') becomes part of the final results. 



4.3 The net correction 5^ and the soft photon exponentiation 

For applications, it is convenient to define a dimensionless radiative correction factor, 

d CTthcoJdEa 



Sa = S{Sa) 
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a = l,m, h, 



(4.38) 
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where d'^a^/d£ is the Born cross section of the process (1.1) and d'^atheor/dS is the theoretical 
prediction for the measured cross section, 



d o"thc 
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-.a \ rfVe d^a-n 
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TT / dt dt 



(4.39) 



One may collect within one expression all the terms, which are explicitly proportional to the 
Born cross section; we call it the factorized part S^^ of the total radiative corrections and 
introduce it in the following way: 
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The net radiative correction factor (4.38) becomes to order 0{a): 
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TT 
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In 5Y-^{£a)) a dilogarithmic term 5inf(£^a) may be found, which controls the multiple soft photon 
emission in the bremsstrahlung process (1.2). An exponentiation of this term leads to a 
significant numerical improvement of the QED predictions [9]: 
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Finally, the radiative correction factor (4.38) with soft photon exponentiation becomes: 
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The factorized correction (4.41) follows from (4.32), (4.34), and (4.35)-(4.37) and is differ- 
ent for the different sets of variables. 
In leptonic variables, it is: 
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(4.46) 
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In mixed variables, we can write (4.43) as follows: 
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In hadronic variables, finally: 
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5 The net radiative correction 5t{ = ^vr + ^r 
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In this section, the net correction will be discussed: 
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(5.1) 



We will use the following abbreviations: 5\cp — ■^i , .^rmx — ^m > 

Before the discussion of the numerical results, some preparations will be performed with 
the infrared finite, hard part of the corrections. This hard part of the cross section has been 
defined in (4.2): 
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From the explicit expressions (3.14)-(3.16) and (4.12) one may see that both terms under the 
integral contain a remnant of the soft photon singularity at Q^ = Ql, while their difference is 
finite by construction. 

The subsequent calculations in this section will be organized in such a way that the in- 
tegrations over Qf or Q\ remain the last ones, respectively. Thus, the compensation of the 
harmless rests of the infrared singularity appears to be quite transparent. In the case of non- 
leptonic variables, in the analytical integrations of the hard bremsstrahlung no special care 
has to be devoted to the soft photon problems. 

In the following subsections, formulae will be presented for the infrared free part of the 
cross section in the different sets of variables. 
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5.1 Leptonic variables 

For the case of leptonic variables, the generalized structure functions Ai{xh,Qh) (i=l,2,3) 
depend on the integration variables yh and Q^ and the structure functions become part of the 
integrand for the remaining twofold integral, which has to be performed numerically. For the 
variables Si = {yi, Qf) and Xi = {y^, Ql), the expression (5.2) is the final result: 



£al 2a'S^ 



dyidQf \i 



Ai{xh, Ql)-^Si{yu Qf, yh, Ql 



dVhdQl J2 

i=l 

J '-'i KUh^i) '- Kyh^i,yh,^h 






(5.3) 



where the radiators Si{yi,Qf, yh,Qh) are given by (3.14) -(3.16). 

The expression for C^^{yi, Qf, y^, Q\) is defined by (4.12) with Q^ = Qf; it is exact in both 
masses m and M: 



C'\yuQlyH,Ql) = T'\yuQlyH,Q 



h) 



?2=Q2 



^737^ + 737^- (5-4) 



The two-dimensional numerical integration in (5.3) has to be performed over the physical 
region (B.34); see figure 35. In the ultra-relativistic approximation, the boundaries are: 



fQl < Ql < min 



'y,S,Qf(l + yi^^ 



(5.5) 



In the numerical integrations, one has to leave out a small region around the phase space points 
with Qf = Ql. There, the integrand is finite but occurs as the difference of two divergent 
terms. 

5.1.1 Discussion 

The radiative corrections in leptonic variables are shown in figures 6-8 for a fixed target 
experiment, HERA, and LEP(8)LHC. At HERA and LEP®LHC, the corrections are very 
similar. For small x, the range of y is reduced due to the condition Q^™™ = 5 GeV^ at a 
given value of x. This tendency is more pronounced for the fixed target experiments. The net 
corrections are smaller there because for the fixed target experiments the S is quite different. 
Further, the muon mass sets the scale in the leading logarithm ln(S'/m^), while at a collider 
it is the electron mass. 

From figure 35 and also from (4.9) it may be seen that for small yi the photon energy 
is strongly bound and only soft photons occur. These corrections are negative. Without 
the soft photon exponentiation, they would even diverge when yi vanishes and at the same 
time xi approaches 1. This may be seen from (4.46) as well. At large yi and small xi 
the factorized correction seems to diverge also; see the squared logarithm in (4.45). This 
behaviour, however, is fictitious. It is compensated by corresponding terms from the hard non- 
factorizing corrections. The steep rise of the corrections at large yi and small xi is completely 
due to the Compton peak, which arises from the small Q\ in the denominator of the photon 
propagator. This is explained in detail in section 8.4. From the Z-exchange diagrams, there 
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Figure 6: Radiative correction Sicp in % for a cross section measurement in terms of leptonic 
variables at a fixed target fip experiment. 
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Figure 7: Radiative correction (5iep in % for a cross section measurement in terms of leptonic 
variables at HERA. 
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Figure 8: Radiative correction 6icp in % for a cross section measurement in terms of leptonic 
variables at LEP®LHC. 

is no contribution to the Compton peak and no steep rise of the corrections at large Q^. This 
may be exphcitly seen in figure 3 of [16]. 

At small values of x;, the corrections are no longer monotonously dependent on xi. In this 
region, the hard photon corrections are dominant. They are not proportional to the Born 
cross section and an influence of the properties of the structure functions at smaller values of 
X may show up. We will come back to this point in section 8. 

The complete QED corrections in leptonic variables have been calculated in the quark 
parton model in [16]. 

5.2 Mixed variables. A second analytical integration 

For the case of mixed variables, it is £m = {Vh, Qf) and the integrations have to be performed 
over the variables Xm = {yi,QD. The cross section (5.2) becomes: 
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(5.6) 



Here, the deflnition (3.12) of the variable Xm is used, Xm = Qf/iyhS). 

In contrast to the case of leptonic variables, in (5.6) one may perform an analytical inte- 
gration over the variable yi. As it is discussed in appendix B.3.3 and may be seen in figure 39, 
the integration region in (5.6) must be split into two regions I and II as long as the condition 
Qf < yhS is fulfilled, which corresponds to x^ < 1- If instead x^ > 1, there is only the 
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region I, see figure 40. Some technical details about the analytical integration over yi may be 
found in appendix D.4. 

In the ultra- relativistic approximation, one gets the following expressions for x^ < 1: 
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(5.12) 



The subtracted part of the cross section also becomes integrated over yi with the aid of 
appendix D.4: 
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The following abbreviations are used: 
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The final formula for (Pa^/dyhdQf can be written in the form: 
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The above expressions again show remnants of the infrared singularity at Qf = Q\. In 
contrast to the case of leptonic variables, one may easily see that the corresponding terms are 
multiplied by kinematical factors, which in the limit Qf = Q\ become explicitly equal to the 
kinematical Born functions Sf, (2.14)-(2.16). This makes the finiteness of the (Pcr^^ obvious. 



5.2.1 Discussion 

The leptonic QED corrections in mixed variables are shown in figures 9-11 for three different 
kinematical situations. As in leptonic variables, at HERA and LEP®LHC they are similar, 
while in the fixed target case less pronounced. 

The cross section depends on Q^ and yh- From figure 39, and also from (4.9), it may be seen 
that for large yh the photon energy is strongly bound; we have only soft photon corrections. 
These corrections are negative there. This is contrary to the case of leptonic variables, where 
this happened at small yi. Without the soft photon exponentiation, they would even diverge 
when yh approaches 1. This behaviour is enhanced, but less than for leptonic variables, if 
Xm approaches 1, too. This may be seen from (4.48). In principle, the Compton peak is 
present in mixed variables. It is considerably suppressed compared to the case of leptonic 
variables, see section 8.4. At small yh, hard photon emission is possible and compensates for 
the negative soft photon corrections. The net correction is positive and remains moderate. If 
in addition x^ approaches 1, we observe that the corrections, if expressed in terms of the S, 
rise steeply. This phenomenon is attributed to the following: the Born cross section does not 
exist for x^ > 1, while hard photon emission does (see also figure 38). So, their ratio explodes 
and the 6 becomes an inappropriate variable. 
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Figure 9: Radiative correction (5mix in % for a cross section measurement in terms of mixed 
variables for fixed target fip scattering. 
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Figure 10: Radiative correction 5mix in % for a cross section measurement in terms of mixed 
variables at HERA. 
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Figure 11: Radiative correction 5mix in % for a cross section measurement in terms of mixed 
variables at LEP®LHC. 
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In sum, the features of the radiative corrections in terms of mixed variables are quite 
different from those in leptonic variables. This is an instructive illustration of the general 
statement that the radiative corrections depend substantial on the choice of variables, in 
which they are determined. In another context, namely for the neutrino-electron charged 
current scattering, a similar observation was made in [60]. 

The complete QED corrections in mixed variables have been calculated in the quark parton 
model in [61]. 

5.3 Hadronic variables. A second analytical integration 

For the cross section in terms of hadronic variables £h = [Vh, Ql), the integrations have to be 
performed over X^ = {yi, Qf) and (5.2) becomes: 
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(5.19) 



In (5.18) we again can perform the analytical integration over yi. As in the case of mixed 
variables, the region of integration in (5.18) must be split into two parts (see appendix B.4.2 
and figure 41). The limits for yi at given values of y^^ Ql, Qf are naturally the same as for the 
mixed variables. The difference is with the Q^, which interchange their roles. 

After the integration over yi, the following expression for d'^a^/dyhQl is obtained: 
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Here the functions 5^' {Qf,yh,Ql 



1, . . . , 3) are exactly the same as those for mixed 



variables [see (5.7)-(5.12)]. The definition of the TC^iQl, Vh, Ql) differs slightly from that of 
the CY^{Qf,yh,Ql)- The integrals which are calculated with the aid of appendix D.4 differ 
also and the result in terms of hadronic variables is: 
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The logarithms L,La,a = l,2,t are defined in (5.15)-(5.16). 

The remarks about the remnants of the infrared singularity and its numerical treatment 
may be simply taken over from the case of the mixed variables. 

The reader may realize that there is one interesting difference between mixed and hadronic 
variables: in the latter case, one may perform also the last integration analytically, thus 
avoiding any numerical part of the calculation. In fact, we have done also this last step, but 
with a different parameterization of the phase space; see section 7.4. 

5.3.1 Discussion 

The leptonic QED corrections in hadronic variables are shown in figures 12-14. They are 
considerably smaller than those which are determined in terms of leptonic variables. From 
figure 41 and also from (4.9), it may be seen that for small yh the photon energy is strongly 
bound; we have only soft photon corrections there. This is the origin of the similarity of the 
gross behaviour of the corrections in hadronic and mixed variables with each other. Without 
the soft photon exponentiation, the corrections would diverge when yh approaches 1. This 
may be seen from (4.50). There is no explicit dependence of the dominant terms on Xh- 
Although, with the factor InQl/m^, one may explain that at higher Xh the corrections are 
slightly more pronounced. Another singularity of the factorizing part of the corrections (4.49)- 
(4.50) is located at yh = 0. As may be seen from the figures, it is compensated for in the 
net corrections. Again, as was mentioned in the case of leptonic variables, this singularity is 
cancelled by a corresponding behaviour of the hard non- factorizing bremsstrahlung. In fact, 
at yh = the corrections vanish even. 
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Figure 12: Radiative correction 6ha.d in % for a cross section measurement in terms of hadronic 
variables for fixed target fip scattering. 
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Figure 13: Radiative correction 5had in % for a cross section measurement in terms of hadronic 
variables at HERA. 
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Figure 14: Radiative correction 5had in % for a cross section measurement in terms of hadronic 
variables at LEP®LHC. 



40 



This may be explicitly seen from the formulae of section 7.4, where a completely integrated 
analytical expression for the corrections in hadronic variables is derived. 

Finally, we should mention that a Compton peak may not be developed since Ql is fixed. 

6 Photoproduction 

The ep scattering process at vanishing photon momentum, Q^ ~ 0, proceeds via the exchange 
of nearly real photons and is called photoproduction. 

The process kinematics usually is described with the Q^ as defined from the leptonic 
variables and the invariant mass W of the compound system consisting of the photon and the 
hadrons, 

Q' = Qf = ih-hf, (6.1) 

W = -iQi+Pi)' 

= M^ + {I - xi)yiS = M^ + yiS ~ Ql (6.2) 

For a fixed value of W'^ << S, the minimal value of Q^ (B.25) may become extremely small. 
With (B.19), W'^ > (M + Tn^)^, it remains non-vanishing: 

Qr\W') - m' ^ / ^ > 0. (6.3) 

The integrated deep inelastic cross section is strictly finite. 

For example, in the reaction ep — > epp at HERA energies the possible values of Q^ extend 
from 10~^^ to 10^ GeV^ [62]. In such a kinematical region, the description of the radiative 
corrections must be exact in both the proton mass M and the electron mass m. For this 
reason, we derive here formulae for the QED corrections in leptonic variables which remain 
valid at extremely small Q^. The chosen integration variables differ slightly from those which 
we used so far: 

£ = £; = iW',Qf), X = T\ = {MlQl). (6.4) 

The invariant mass M| of the hadronic system is 

Ml = -{Vi + Qhf = M' + ynS - Ql (6.5) 

It will be convenient to treat the infrared singularity in a lorentz system where the above 
mentioned compound system is at rest: p2 + ^ = 0. The necessary kinematical relations are 
derived in the appendices A. 2 and B.2. In the following it will be described how the three 
different contributions to the corrections are treated: the infrared divergent correction ^g^^, 
the factorized hard part ^^^rd) ^^^ the finite rest of the hard bremsstrahlung contribution. 

The removal of the infrared divergence proceeds technically as it is described in section 4.2 
but will be performed in the rest system introduced here. For the soft part of 5^^ one gets 
again the expression (4.26). The velocities /5j, which are introduced in (4.20) are now to be 
expressed by (A. 7) with the invariants (A. 6). The result is: 
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4m2M2 

41 



4m2M2 , ^ 

(6.6) 



212 ■ 



Further, for 5$ the exact definition (D.22) has to be used in terms of the momenta (A. 7). 

Integrating (4.13) threefold as it is explained in appendix D.3.2 one obtains the following 
expression: 



CrdiVl^Ql^) 



-In h m ; 



2P 



L, 



(6.7) 



where 



Ln = In 



(3 + 1 



(3 = 



Am? 



(6.8) 



Further, we need the exact QED vertex correction, which contains ^vert [6]: 

5vcrt(/3) = 



-(-■^-'■'^(^■^-04--^ 



2/9 



Ijnhl-::^ + Li2 -: - Li2 



-/3J 



(3^ 



1-/5. 



l + (3 



(6.9) 



Finally one has to take into account also the contribution from the anomalous magnetic 
moment of the electron [63]. This vertex correction is needed only for the virtual photon 
exchange, i.e. one may neglect the axial vector part in the matrix element (2.1): 



7m -^ SvcTtiPhi, + imVi,nom{P) 



{h + /C2)^ 



'^anoml.PJ 



hi 



(6.10) 
(6.11) 



The cross section contribution from the anomalous magnetic moment is" 



(Pa^ 



2a^S m^ 



dyidQ^i XsQf xiQt 



Vanom(/3) \2(3^yfxiF,{xi, Qf) - (2 - mfF,{xi, Q\ 



(6.12) 



Combining the three factorizing corrections, one gets for the factorized part 5vr the fol- 
lowing expression: 

SyRiVh Ql ) = <^vcrt + (^hard + ^soft 

X ( n2\^ 1 1 l+/5i ^ 1 1 1 + /52 ^ <, 



'$ 



2' ^ 2(3 



T 1 4/9' , T • /l+/5^ 
Lflln-::: h LI2 



-/3i 



P' 



1-/9, 



Li-; 



• 1-/9 ^ 
1 + /9, 



where 



A ( n2\ 91 W^-{M + m^f ( l+(3\ 
Omiivi, Ql) = 2 In j== I ^^ L/3 - 1 I . 



m 



y/W 



2(5 



(6.13) 



(6.14) 



^The corresponding equation (47) in [6] was not correct. 
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Finally, the infrared free part of the cross section (1.1) can be written as follows: 



d^al 2a^S 



dyidQf Xi 



i=l 



1 

'Qi' 



M^h, QD^r^Siiyu Ql Vh, Ql) 






(6.15) 



where Si{yu Qf, Vh, Ql) {i = 1, 2, 3) are given by (3.14)-(3.16), and C^^{yi, Qf, yu, Ql) is defined 
by (5.4). We would like to remind that these expressions are exact in both masses m and M 
for the photon exchange contribution. 
In sum, the cross section is: 



(i^o-R rf o-anom , d^a^ , rfVe 



-)2 ~ j„ j/n2 ' j„ j/n2 "'" j„ j/n2 l ^-^P 



dyidQf dyidQf dyidQf dyidQ 



-S-miiS) 



1 + ^ 

TT 



SMS)-6i^f{S) . (6.16) 



Here, we should also mention that the structure functions must have a certain behaviour in 
the limit of small Q^. From the relation between the total cross section of photon-proton 
scattering and the structure function F2 at small Q^ [62], 

^1{W) = ^F,{x,Q')\q.=o, (6.17) 

one concludes that the structure functions should vanish with vanishing Q"^. One possibility 
is to multiply them by a global suppression factor [64]: 

J^5^(a;,Q2) ^ [1 _ exp{-aQ^)]J^^^^{x,Q^), a = 3.37 GeV-^. (6.18) 

In our numerical results, we follow this prescription. 

6.1 Discussion 

The QED corrections in the case of photoproduction are shown in figure 15 for HERA energies. 
The structure functions [54, 65] are taken from [55] and [66]. 

As was already discussed in the section on leptonic variables, soft photon emission is located 
near yi = 0. With rising yi more and more hard photons may be emitted and the influence 
of the structure functions via the non- factorizing hard corrections becomes more pronounced. 
It is further interesting to observe that at the extreme small values of xi as discussed here 
the radiative corrections vanish at yi = 0. For the anomalous vertex correction (6.12) and the 
exponentiated soft photon correction (6.14) this follows immediately from the limits in (6.8): 
/5 — > cx) and Lp -^ 2//9. Here, it is important to notice that the Qf becomes small even 
compared to the electron mass. In fact, with the aid of (6.2), (6.3) may be well approximated by 
Q2in"i ^ m'^yf << m^ for small yi. Further, the f^^^rd of (6-7) vanishes explicitly. This reflects 
the general statement that for sufficiently small yi energetic photons cannot be emitted. For 
the rest of the corrections, we see no trivial argument why they should vanish or compensate 
each other; but they do. A further interesting feature is of kinematical origin. At very small 
xi, the Qf is also bound strongly. As may be seen in figure 35, then the allowed values of yu 
are bound from above. In the ultra-relativistic approximation this property gets lost. 
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Figure 15: Radiative correction for photoproduction at HERA with two different sets of struc- 
ture functions. 

7 An alternative treatment of the phase space 

So far, one unique phase space parameterization has been used for all three different sets of 
external variables. The infrared problem was treated in the proton rest system. Compared 
to this, for leptonic variables the rest system which was introduced for the treatment of the 
photoproduction process lead to considerable simplifications. 

Here, we change to a completely different phase space parametrization, which is specially 
adapted to the case of hadronic and Jaquet-Blondel variables. The leptonic degrees of freedom 
will be integrated over at the very beginning of the calculation. 

7.1 Hadronic variables. The phase space 

Wherever possible, the ultra-relativistic kinematics will be used in this section. We start from 
the slightly rewritten expression (2.38): 



dk2 dk 
2k2 2A;0 



S\ki +pi-k2-p2- k)d%6{pl + Ml)dMl 



(7.1) 



Now, the integration variables we are interested in will be introduced. For that purpose, we 
consider the photon and the final state electron to be a compound system. The corresponding 
rest frame R is defined by the three-momentum relation k2+k^ = 0. The integration variables 
k2, k may be replaced by the following variables: 



the invariant mass of the (76) compound system r = —{k2 + k) 



2. 



44 



• the photonic angles 'dpi., fR in the rest frame R of this compound system. 
The phase space of the variables k2, k may be parameterized as follows: 
dk2 dk 

Here, the phase space element of the (7e)-compound system is introduced: 



/^^5'(A:i+Pi-A;2-p2-A:) = j dr5{T + K^)dT,,. (7.2) 



In the rest frame R of the compound system, (7.3) may be rewritten as it is explained in 
appendix A. 4. 3: 

7A(r,m^0) ri r2. 

dT^e = z / dcosi3R / dipR, 

8t J-1 Jo 

^A(r,m2,0) = T-m'^. (7.4) 

In a next step, the hadronic momentum transferred Ql is introduced into (7.1), and the 
corresponding (5-function under the p2 integral is exploited. In appendix A. 4. 4, it is shown: 

/ d'p2S{pl + Ml)5 Uk, + pi - p2f + r] 5 [qI - {p2 - pif] = -^. (7.5) 



2VAi 

Then, with the relation yhS = Ql + Ml — M^ from (6.5) one gets the identity dQ\dMl = 
SdQidyh- This introduces the last of the hadronic variables into the integration measure: 

-TT-^ S* /" f It 777 f 

r = / dyhdQl / dr- / d cos i^RdcpR. (7.6) 

4v Aq J J 4:71 T J 



The physical region of {yh, Ql, t) is derived in appendix B.4.4. 

7.2 A twofold angular integration 

In hadronic variables, the double differential cross section of the process (1.2) becomes: 



dyhdQl SQl . .^^ 



j dTY^A,{xh,Ql)S,{yh,Qlr), (7.7) 



where 

Si{yh,Ql,r) = -^^ J d cos i^RdipRSi{yh,Ql,T, cos i^R,ipR). (7.8) 

The functions Si{yh, Ql, t, cos'^r, ipr) arise from the functions Si{S,X, z\, z^) after the follow- 
ing substitutions for t/^, Q^, and z\ in (2.40)-(2.42) in accordance with appendix A. 3 and with 
the definitions (2.34), (2.43): 

Ql = QI + Z2-Z,, (7.9) 

zi = ^(Ql + T + m^-^icos^R^, (7.10) 

I / ■TOO \ ^ I 

1-yi = -(l-|/ft)fn \ + -^y\ricosi9pCosi^R + smi^pSmi^RsmipR), (7.11) 
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where 



cos-i^r- 



S{Ql - z,) - ynS{Ql + r + m^) 



A^yAi 



(7.12) 



After performing the above insertions, the Si may be integrated over the two photon angles 
with the aid of the tables of integrals of appendix E.l. In the ultra-relativistic approximation, 
the result is: 



Q 



1 



Z2 



(L. 



Si{yh,Ql,T) 

S2{yk,Qlr) = S'{2{l-y,) 



S3{yh,Ql,r) 



where 



1 r — m 
-(L.-2)- 

Z2 



+ T + - I 1 - ^ I W, 



1 1 

4 + 2 
It — rn^ 



Ql 



+ ^[l-(l-l//.)(W-3)] 



s\2Ql{2-y^) 
Ql 



-(W-2) 

Z2 



4 r2 
1 

It — irP' 



or 

Qh 

|(L.-3) 



(7.13) 
(7.14) 



{2-yj,)U--{l-yh){7-yh) 



+ 



Q'r 



2-{U-2)l2-y,,+ 



4 r2 



yft L 



(7.15) 



U 



In 



m'^T 



(7.16) 



g^ = Q2+r-m2. 



(7.17) 



The infrared singularity is located at 2:2 = and has to be treated appropriately. 



7.3 The infrared divergence 

As was done in section 4.1, the following ansatz will be used in order to separate the infrared 
divergent parts in (7.7): 



dyhdQl 



-^ E A{x,, Ql) \ Sf{y,, Ql)6'^{y,, Q 



h i=l 



+ dT 



S,{y,, Ql t) - Sf{yn, QDJ^'^iyn, Ql r 



(7.18) 



Here, we define the S^{yh,Ql) in complete analogy with (4.14) as the sum of an infrared 
divergent soft and a finite hard contribution: 



d^k 



2 f a k 
SlUyH,Qle) = -j^,:F'\Qlz,,Z2)e{e-h: 



Sh.rdiyh,Qle) 



4 v^ rdV 



(7.19) 



TT^ 



S 



j^:F'\Qlz,,z,)e{k'-e) 



dT / dcos^Rd^Rj^^\Q^ , zi, Z2)e{k^ - t). (7.20) 

Atx t 
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The J-'^^{yh, Q\, t) will be determined below. The expression for (5g^t(|//i, Q\, e) is the equiva- 
lent of (4.10) and Si^rAVhiQ'hi^) follows immediately from (4.14). Further, 



^ (Qft,^l,^2) 



Z1Z2 



m 



1 1 

^^1 ^2y 



(7.21) 



The hard part of the infrared divergent correction will be integrated first. The integrals over 
the photon angles in (5hard(l/^? Qhy ^) ^^Y be performed with the aid of (7.10) and of the table 
of integrals in appendix E.l: 



S^ZdiVh, Ql e) 



47r T J 
1 11 

- (W - 2) - —K + -, 

Z2 Qi r 



(7.22) 



(7.23) 



and L,- and Q^ are defined in (7.16) and (7.17). 

The boundaries for the integration over r may be found in (B.94). The potentially infrared 
divergent term in S^j.^{yh, Qh, e) contains a dependence on 1/2:2 and is integrated with the aid 
of the table of integrals E.2: 



r <*- ,L, - 

Z2 Z2 
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- 2Li2 1 - - 


\Ql(i A] 

_2me \yh J 

I \ 

- -Li2(l). 
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/I 

\yh 


-1 



Vh, 



(7.24) 



With the aid of appendix E.3, the exphcitly finite part of 5^-^^{iih,Qle) may finally also be 
integrated. We will not quote the result separately, but rather treat this contribution together 
with the finite hard cross section in the next section. 

Now, the infrared divergence in 5^^^^^{iih,Q\,() will be treated. For this purpose, the R 
system of the foregoing section is used. As derived in appendix A. 3, the photon energy 
becomes: 
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,0,R 
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R\ 



T — m 



Z2 



< e. 



2v^ 2^ 
Besides A;-^, we will need the four momenta /cf , /cf^ in the limit Z2 — * 0, see A. 9: 



(7.25) 
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= 1^^ ~ Z2 - 


- 0, 


A;?'« - 


_ Ql 
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, Ql 
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4m^ 



(7.26) 



(g^ + 2my 

Since the final state lepton is practically at rest there, the following isotropic relation holds: 

-2A;A;2 ^ 2mA;°-^. (7.27) 

It is this property, which simplifies the angular integration. In (4.15), it is not necessary to 
introduce the Feynman parameter which finally lead to the complicated function 5*$ of (D.22). 
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Instead of (4.16), one gets: 



2(27r) 
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where one has to use Q^ -^ Q\. Instead of (4.23) and (4.24), it is now 
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(7.28) 



(7.29) 



(7.30) 



The integrations in (7.29) may be performed with the aid of appendix D.2: 



C,,{y,,Qle) 



>IR 



2e 



P^" + In — 
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L^ + Lh - Li2(l) 



1, 



(7.3i; 



where V^^ is introduced in (4.25), and 



Lh = In 



Ql 



m,^ 



(7.32) 



The cut-off parameter e may be seen to be cancelled against the corresponding terms in (7.24). 
In order to get an infrared finite and well-defined cross section the photonic vertex correc- 
tion (5vert has to be added to the radiative cross section (7.18). It is given by (4.34), again with 
the replacement Q^ -^ Ql- In the sum, the infrared divergence cancels: 



CtiVh, Ql e) + 5^UQl) = 2 (Lh - 1) In ^ - ^L^ + ^L^ - 1. 



(7.33) 



7.4 A third analytical integration 



In (7.18), the last remaining integral is that over r or, equivalently, over Z2. Since the structure 
functions depend on hadronic variables which are here at the same time the external variables 
the last integral may be performed analytically. For this purpose one has to integrate the 
functions Si{yh,Qh,T) as given by (7.13)-(7.15). The cross section may be reordered such 
that the infrared singular terms in the finite hard cross section explicitly compensate each 
other: 
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(7.34) 
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The resulting cross section is rather compact and the last integrations over r may be easily 
performed with the aid of appendix E.3. 
The net cross section is: 



(Pa (faB 



dyhdQi dyhdQ 



exp 



-Sini{yh,Ql) 



+ ^ E 7^A(x., QDS^ivH, Ql). (7.35) 



The (5inf(l//i, QD and the Si{yh,Qh) are defined in (1.10) and (1.11). 

We would like to remark that these purely analytical formulae for the leptonic QED cor- 
rections in hadronic variables to neutral current deep inelastic scattering are an extremely 
esthetic result. They are published here for the first time. Numerically, the corrections agree 
completely with those derived within the approach of section 5.3. 

As was mentioned above, the corrections vanish aX y ^ 0. That this is really the case may 
be seen by inspecting the analytical expressions. 

7.5 Jaquet-Blondel variables 

From hadronic variables, one may easily change to Jaquet-Blondel variables (1.8) in the phase 
space (7.6). The expression for Qjb in terms of hadronic variables is derived in appendix B.5: 

1/jB = Vh, (7.36) 

Q% = Ql-^{r-m\ (7.37) 

The physical boundaries in the phase space integral (7.6), 

r = — / dy:iBdQ% / dr- / dcos'^ndipR, (7.38) 

are also derived there. From the above relations, it may be seen that the Jaquet-Blondel 
variables and the hadronic variables are related in a one-to-one correspondence, with r being 
a parameter. 

In Jaquet-Blondel variables, the double differential cross section of process (1.2) reads: 



(i^o-R 2a'^ 



, ,n^ o f drJ2Mx,,Ql)^S,iyjB,QJB,r). (7.39) 

The iSj(|/jB, QjB) ''") iiiay be trivially obtained from (7.13)-(7.15) with the relations 

Vh = yjB, (7.40) 

Ql = QJB + -r^^ir-m'). (7.41) 

1 -yjB 



The explicit expressions are [46]^: 
5i((5jB,l/JB,r) = QJb 



-(L.-2) + f^' 
Z2 4r^ 



1 - SyjB . (^ ^2 yjB \ .„ ,. 



^Thc equation (7.44) corrects the equation (17) in [46]. The error was created in the manuscript. Fortran 
program and numerical results were not influenced by this. 
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(7.43) 
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•. (7.44) 



The integral over r has to be performed numerically as long as one cannot neglect the 
difference between Qjb and Ql - the structure functions are depending on Ql and thus on r. 

The integrand in (7.39) is infrared singular at r = m^, which corresponds to Z2 = 0. The 
treatment of the infrared singularity may be performed in full analogy with the foregoing 
section. In the infrared singular point, Z2 -^ 0, the hadronic and the Jaquet-Blondel variables 
agree and the 5^^^^ in these variables do so: 
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The (^hardd/JB, QjB) ^) differs slightly from the hadronic correction: 

4 v^ rdV 



'^hard(l/JB,<5jB)e) 



It — irP' 



dr- 



An T 



dcosdRd^Rj^^^{Q%,zi,Z2)e{k^ - e). (7.46) 



In fact, it is 
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(7.47) 



Using the results of the calculation of the hadronic correction, it is sufficient to explicitly 
integrate the difference term in (7.47) over ip^ and cos'dji with the aid of appendix E.l. Then, 



Cd(l/JB,gL,e) = J dTj^'%j^,Qj^,T)e{k' 
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The Lt- and Q^ are defined in (7.16) and (7.17), respectively; furthemore the substitutions 
(7.40) and (7.41) have to be used. The integral over r may be performed with appendix E.2 
for the potentially infrared singular terms and with appendix E.3 for the others; but now 
within the boundaries (B.106). 
The result is: 
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The photonic vertex correction (5vcrt is given by (4.34), but with the replacement Q^ 
The net factorized part is 
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The finite part of the radiative cross section in Jaquet-Blondel variables is 
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The net cross section is: 
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7.5.1 Discussion 

The leptonic QED corrections in Jaquet-Blondel variables are sliown in figures 16-18. Tliey 
behave quite similar to the corrections in hadronic variables. At small y, this follows from the 
fact that the variables become equal there. Thus, the vanishing of the corrections at y^ — » is 
understood from the hadronic case. As was discussed with the mixed and hadronic variables, 
the soft photon corner is at |/jb approaching 1. From (7.55) one may see that the dominant 
soft photon corrections are negative and the more pronounced, the larger the xjb. This makes 
some difference to the hadronic variables. Again as in all the previous cases, the corrections 
for fixed target scattering are different from the collider cases due to the smaller value of S 
and to the difference of the muon and electron masses. The collider results are quite similar 
to each other. 



8 Discussion 

In the foregoing sections, we discussed some of the basic features of the leptonic QED correc- 
tions to the neutral current deep inelastic ep scattering in terms of different sets of kinematical 
variables. 

Basic features of the magnitudes of the corrections and their behaviour as functions of the 
kinematics have been discussed there. 



-10 7 



-20 



-30 



-40 7 




Fixed target /zp 



-50 ^ 

0.0 0.2 0.4 0.6 0.8 1.0 

y 



Figure 16: Radiative cross section for fixed target deep inelastic neutral current scattering in 
terms of Jaquet-Blondel variables. 
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Figure 17: Radiative cross section for deep inelastic neutral current scattering at HERA in 
terms of Jaquet-Blondel variables. 
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Figure 18: Radiative cross section for deep inelastic neutral current scattering at LEP^LHC 
in terms of Jaquet-Blondel variables. 
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In this section, we discuss some additional features of the corrections: 

• the dependence of the corrections on the structure functions; 

• comparison with the leading logarithmic approximation; 

• photonic cuts in leptonic variables; 

• the Compton peak; 

• other radiative corrections. 

8.1 The influence of structure functions on the QED corrections 

In figures 19 and 20, the radiative corrections in leptonic variables at HERA and LEP®LHC 
are parameterized with two different sets of structure functions [54, 55, 67, 68]. Both sets 
were derived before HERA started operation. If there are deviations between predictions with 
different structure functions, then they are due to the contribution from hard photon emission. 
For soft photons, the correction factorizes and the structure function dependence drops out in 
the corrections, when expressed in form of the ratios 6. In experimentally unexplored regions 
the predicted structure functions rely heavily on the predictive power of perturbative QCD. 
Especially for small values of x problems may arise. These are reflected in rising deviations of 
corrections with different structure functions when x becomes smaller. While, at intermediate 
and larger values of x, the dependence of the corrections on the choice of structure functions 
becomes minor. For this reason we changed the presentation of the corrections and chose the 
set of variables Q'^,x, thereby losing sensitivity to y. 

Figures 21 and 22 show the same corrections as figures 19 and 20, respectively. But now 
with structure functions [69, 70, 67, 71] whose derivations took into account HERA data at 
low X. The disagreement of the predictions at small x remains but is much smaller. 

A similar comparison for the other sets of variables is shown in figures 23-25. The structure 
functions [54, 55, 67, 68] which were not yet improved with the HERA data are used. The 
dependence of the corrections in mixed variables on the choice of structure functions is quite 
substantial but less pronounced than in the corresponding case in leptonic variables. For the 
hadronic and Jaquet-Blondel variables the influence is reasonably small. This follows from the 
smaller contribution of hard photon emission. The net corrections are negative. This proves 
the dominance of the factorizing soft and virtual photonic contributions. Their dependence 
on the structure functions cancels when 6 is calculated. 

The dependence of the QED corrections on the structure functions is not only influenced 
by different estimates of the QCD predictions for the parton distribution evolutions. The 
corrections show additional dependences on the low-Q^ behaviour of the structure functions, 
the choice of Q^ in the argument of the parton distributions, soft photon exponentiation, and 
the running of the QED coupling constant. 

We use the corrections in mixed variables for an illustration of this. The Fortran program 
TERAD91 [41] allows us to switch on and off the above indicated dependencies of the cross 
sections with the following flags: 

(i) For IVAR=0, the low-Q^ behaviour of the structure functions remains untouched; for 
IVAR=2, the modification (6.18) is applied. 
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Figure 19: Comparison of the radiative corrections (5iep(a;, Q^) at HERA in terms of two dif- 
ferent sets of structure functions. 
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Figure 20: Comparison of the radiative corrections (5iep(x, Q^) at LEP®LHC in terms of two 
different sets of structure functions. 
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Figure 21: Comparison of the radiative corrections (5iep(a;, Q^) at HERA in terms of two dif- 
ferent sets of HERA-improved structure functions. 
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Figure 22: Comparison of the radiative corrections (5iep(x, Q^) at LEP®LHC in terms of two 
different sets of HERA-improved structure functions. 



56 



20 



:5™.x[^"] 



-20 




HERA ep 



MRS D'_(DIS) 
CTEQ2D (DIS) 



-40 



1 10 10 ' 10 



10 * 10 

Q'[GeV'] 



Figure 23: Comparison of the radiative corrections S„^ix{x, Q"^) at HERA in terms of two 
different sets of structure functions. 
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Figure 24: Radiative corrections 6had{x, Q"^) at HERA in terms of two different sets of structure 
functions. 
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Figure 25: Radiative corrections 6jb{x, Q"^) at HERA in terms of two different sets of structure 
functions. 

(ii) For ITERAD=0, the structure functions are assumed to depend on Qf. A correct choice 
is ITERAD=1, which makes the structure functions dependent on Q\. 

(iii) With the settings IEXP=0,1 the soft photon exponentiation as it is introduced in (4.43) 
is switched off or on. 

(iv) For IVP0L=0,1 the running of the QED couphng constant with the Q^ of the t-channel 
momentum flow is switched off or on. 

In the figures, the flag settings are quoted as follows: (lVAR,ITERAD,IEXP,IVPOL). 

Figure 26 shows the dependence of the corrections on both the running of the QED cou- 
pling constant and the soft photon exponentiation. Both corrections are small and tend to 
compensate each other for large values of y. For Xm = 0.5 and the maximal value of Q^, the 
values are: -34.4% for (0000), -29.8% for (0010), and -32.7% for (0011). 

Figure 27 shows the dependence of the corrections on the choice of the Q^ in the structure 
functions under the integral for the hard photon emission. There is a strong dependence 
when hard photon emission contributes substantially to the corrections. We see also a faking 
behaviour of the corrections at smaller values of Q^, which is more pronounced at larger values 
of X. At Q^ = 5 GeV^, the parton distributions are artificially frozen since below this value 
the evolution may not be controlled. Thus, there is no room for a variation of Q^ at all, and 
both flag settings become equal. 

Finally, figure 28 shows the most realistic choice of the flags mentioned. As is evident from 
the above discussion and a comparison with the foregoing figures, the low-Q^ modifications 
are extremely important for a correct description of the corrections at small Q^ and large x. 



58 



100 



75 



x = 0.5 



50 - 



25 



-25 



-50 



HERA ep 



.x[%] 



10^*1 10"^! 10"^ 1 10"^ 




11 



10 10 ' 10 ' 10 ' 10 ' 10 



Q [GeV^] 



Figure 26: Comparison of the radiative corrections Sniix{x,Q'^) at HERA with two parameteri- 
zations of structure functions. Flag settings (^IVAR,ITERAD,IEXP,IVPOL^ are explained in the 
text. 
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Figure 27: Comparison of the radiative corrections 5mix(a;, Q"^) at HERA with two parameter- 
izations of structure functions; the flag settings are explained in the text. 
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Figure 28: Comparison of the radiative corrections S^-aix{x, Q^) at HERA with two parameter- 
izations of structure functions. Flag settings are explained in the text. 

To summarize this part of the discussion, at least in certain regions of the phase space it 
is extremely important to ensure an iterative improvement of the various components of 
the description of deep inelastic scattering, including the interplay of structure functions and 
QED corrections. 



8.2 A comparison with leading logarithmic approximations 

An interesting and practically important question concerns the accuracy of the leading logar- 
ithmic approximations (LLA). For a comparison, we use the Fortran program HELIOS [36] 
with the same structure functions as in TERAD. For the special purpose here soft photon 
exponentiation, other higher order corrections, the low-Q^ modifications etc. are excluded. 

The comparison is performed for leptonic, mixed, and Jaquet-Blondel variables and shown 
in figures 29 - 31. For the comparison, the structure functions [67, 72] and the following flags 
in TERAD have been chosen: (IVAR,ITERAD,IEXP,IVPOL) = (0100). 

Over a wide kinematical range, the leading logarithmic approximation works quite well 
and is completely sufficient for a description of the experimental data. For intermediate x in 
leptonic and small x in mixed and hadronic variables, the accuracy of the LLA is best. 

One may also conclude from the figures that a complete 0{a) calculation will become 
necessary if the experimental accuracy reaches the level of one per cent or better. 
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Figure 29: Comparison of radiative corrections S\cp{x,Q^) at HERA from a complete 0{a) 
and a LLA calculation. 
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Figure 30: Comparison of radiative corrections Smixix, Q"^) at HERA from a complete 0{a) 
and a LLA calculation. 
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Figure 31: Comparison of radiative corrections Sjb{x,Q^) at HERA from a complete 0{a) 
and a LLA calculation. 

8.3 Photonic cuts 

An instructive illustration of some basic features of the QED corrections is figure 32^. It shows 
the dependence of the correction 6\cp on yi in presence of cuts on the photon momentum which 
are introduced in appendix B.1.2. 

To be definite, we chose x = 0.001 and HERA kinematics. The quark distributions are [65, 
66] . It is assumed that photons may be observed if they have an energy larger than Em = 1 
GeV and are inside a cone with an opening angle of ^™^^ = 0.1 rad. The cut conditions are 
inactive as long as the photon energy is smaller than Em- The left-hand side of the second 
of (B.15) corresponds to this. As long as z™"^ is negative, the cut is fictitious. It starts to be 
influential if, in the ultra-relativistic limit. 



, Em 

yi = yh + -^■ 
Since yh is integrated over, this condition weakens to the following one: 



U] 



yi 



Em 



0.033. 



f8.2) 



At this value of yi the rise of the cross section is sharply interrupted. The collinear initial and 
final state photons may get an energy larger than Em- If so they become cut when being 



^Thc figure was made with data from TERAD91. Later, it was reproduced by the Monte-Carlo program 
HERACLES [32, 38]. We should like to thank H. Spiesberger for his careful analysis of the figure, which we used 
when preparing this section. 
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Figure 32: Radiative cross section for deep inelastic neutral current scattering at HERA in 
leptonic variables with photonic cuts; x =0.001. 

emitted inside the cut cone. The colhnear final state photons escape from the cut cone and 
become non-observable if the electron scattering angle 6'^, 
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becomes larger than ^™^^. The value of yi above which this is the case and the collinear final 
state photons become more and more unconstrained by our cuts is: 
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'l-cos9'jEe 
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0.084. 
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Above this value, the contribution from final state radiation to the cross section rises. 

Further, the emission of photons due to the Compton mechanism grows. This contribution 
peaks at large yi and small xi (see the next section). It is due to the t-channel exchange 
of photons with practically vanishing Q^. Thus, it comes mainly from events with real pho- 
tons whose transverse momenta are balanced with those of the scattered electrons and the 
longitudinal ones are opposite. If now yi becomes larger than the value 



yi = i 



)Ee 



(l-cos^^ 

(1 -COS^^)Ee + (l + COS ^; 



_ — _ = 1-0.084 = 0.916, 
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these photons enter into the cut cone and will not contribute to the measured cross section. 
The Compton peak starts to become cut. 
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8.4 The Compton peak 

The Compton peak has been discussed first in [2]. The corresponding cross section enhance- 
ment is due to energetic final state radiation from the lepton. It is related to large yi and 
small Xi. As will be shown here, the origin of the Compton peak is deeply connected with the 
lower bound on the Q\ which is reached when the integral over Q\ is performed in order to 
get a cross section in leptonic variables. The extremely small Q^ of the virtual photon makes 
the kinematical situation resembling the scattering of light at matter. 

As may be also seen from the figures, the Compton peak occurs in leptonic and mixed 
variables, while it is absent in hadronic and Jaquet-Blondel variables. The reason is that in 
the last two sets of variables, the Q\ is an external variable and thus not free to vary down to 
the kinematical limit. 

The Compton peak arises from the integration of the photon propagator dQ\/Q\ over Q\ 
in the contribution from the 7 exchange in formulae (5.3) and (5.6). We will indicate here the 
essential steps only. For the structure functions, Bjorken scaling will be assumed: 



F,{xh,Ql) - F,{xh). 



\.7) 



Further, for an integration over Q\ at fixed Xh the following change of integration variables is 
necessary: 



di/hdQl = dxhdQ 



01 

^xlS- 



The variables Xh and Q\ vary within limits which are derived in appendix B.2.5. 

Using (8.7) and (8.8), the 7 exchange part of the radiative cross section may be rewritten 
as follows: 
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(8.10) 



The correction (5vr may be found in [5, 16]. For a discussion of the Compton peak, only 
the first two terms under the integral in (8.10) are relevant. Examining their variations 
as functions of Qf one may conclude that they are of two types. From the integrands 
(1; ll\fC[\ \ I \fC'^dQ\l Qf, a ln((5^™"'/(5^™'''') emerges, while from the integrands (1; \l\fC[] 
\ j \/G2)dQfl Q\ a l/Ql™"^. The latter, however, are screened by factors rr? or M^. 

In [5, 16], the following formula has been derived after integration of (8.10) over Qf in the 
ultra- relativistic approximation: 
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Y ml J Ji z — \ 

dz— [f/e(7; 1, i-yi) + Ueir, yi - i, -i)] F^izxi) , (8.11) 
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From the exact boundaries of appendix B.2.5 one gets the following relations: 

Qr-'iyuQlxn) ^ ^Ql (8.14) 

Xl 

Qr\yuQlxn) - .^^ ''f'f ^ M\ (8.15) 

Therefore, one of the logarithmically peaking integrands in (8.10) is: 
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This contribution is explicitly seen in (8.12). Two additional logarithmically enhanced con- 
tributions may be found in the last term in (8.12). The above mentioned screening of the 
other potential sources of kinematical peaks in the variable Q\ may also be seen in (8.11) - 
the inverse powers of M^ are absent. 

The argument of the logarithm at the right hand side of (8.16) may become very large 
at xj ~ 0,1// ~ 1 {xh > Xl). This is the origin of the Compton peak. We would like to 
stress that the cross section depends on the proton mass M on purely kinematical grounds. It 
has nothing to do with assumptions in the quark parton model. In the above considerations 
we used one implicit assumption. Namely, for the argument to hold, the parton distribution 
must be a slowly varying, non- vanishing function at smallest Q\. This is exactly, what they 
aren't. We remember here about the discussion of (6.3). Thus, the Compton peak is at 
least partly regularized by this damping of the quark distributions. In view of lack of a 
reliable model for the parton distributions in the region of smallest Ql it is not evident which 
theoretical prediction of the radiative corrections is best. Maybe one should look inversely 
onto the problem: a clean measurement of the Coulomb peak could teach us something about 
the parton distributions which give rise to it [73]. We do not think that there are strong 
arguments to favor Aqcd as a scale of damping of the parton distributions at small Ql as was 
proposed in [74]. 

So far, we discussed the case of leptonic variables. In section 5.2.1, it was mentioned that 
the Compton peak in mixed variables is much less pronounced but present. The cross section 
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in mixed variables has been integrated once, over yi, in (5.7)-(5.17). Only in one of the two 
regions, in region I, the Q\ may become small. There, again the integration starts from Q\ 
The minimal value of Q\ has been determined in (B.62). With the aid of (B.22), in the 
ultra-relativistic limit in mg it is yf^^^iQf) = 1 — M'^Qf/S'^ and 

Qr"" = QUi-yr^' + yk) 

= XraVlil + rXm)S ^ XmVlS. (8.17) 

The corresponding kinematical singularity is located at x^ —* 0,yh ^ 0. The extreme value 
of Ql for mixed variables is scaled by a factor of the order of S/M"^ compared to the case 
of leptonic variables and the peak is much less pronounced in the experimentally accessible 
kinematical ranges. 

8.5 Summary and outlook 

In this article, we gave a comprehensive presentation of the 0{a) leptonic QED corrections 
to deep inelastic ep scattering. The detailed analysis of the kinematics of deep inelastic 
scattering is of a general interest. It has been performed exactly in both the electron and 
proton masses in order to scope both with extreme kinematical situations and the characteristic 
singularities of the totally differential cross section. Then, for several sets of kinematical 
variables semi-analytical integrations have been performed. The remaining two integrals for 
leptonic variables have to be taken numerically while for mixed and Jaquet-Blondel variables 
one integral remains. The hadronic case has been solved completely analytically. Most of the 
results of this article have not been published before. 

It was not the intention to cover all relevant radiative corrections. 

Here we would like to mention which corrections have been left out of the presentation: 

• Electroweak one loop corrections [16]. These have to be treated in the quark parton 
model. 

• The QED corrections which are not related to the lepton legs - the lepton-quark inter- 
ference and the quarkonic radiation. Their calculation rests also on the quark parton 
model. The corrections are known in leptonic variables [16] and in mixed variables [61]. 
In other variables they are under study. 

• Within the LLA approach, there exist predictions for QED corrections also in other 
variables; for an overview, see [78]. A semi-analytical, complete 0{a) calculation in 
these variables seems to be nontrivial. 

• A discussion of the Monte-Carlo approach has been left out completely. For this topic, 
we refer to the overviews given in [23, 24]. 

• Further, there is the wide field of higher order corrections. Their importance depends 
on the variables chosen. For leptonic variables. They have certainly to be included in 
certain regions of the phase space if an accuracy of one per cent or better is aimed at. 
The soft photon exponentiation is part of this and the hard photon contributions beyond 
0{a) are certainly needed in LLA. Partial solutions are obtained in [37, 76, 77, 78]. 

• Aiming at an accuracy of 1% or better, one has to concern about QCD corrections to 
the parton distributions. This has been done in the HECTOR project [78]. 
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• Finally, one should mention that for many physical applications charged current scat- 
tering has to be described with the same accuracy as the neutral current one. A lot of 
theoretical work remains to be done here. So far we may refer only to [18, 32]. 

To summarize, in this article we treated the numerically largest QED radiative corrections 
semi-analytically in several of the most important kinematical variables. This gives an im- 
portant tool for a satisfactory description of fixed target and HERA data and the possibility 
to estimate the corrections at much higher beam energies. The most important next (or 
alternative) steps have been indicated. 

Note added in March 1996. 
We corrected misprints in the following formulae: (2.25), (A. 23), (A. 26). 
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A Kinematics and phase space 

A.l Kinematics and phase space for section 3 

Throughout this section, all definitions and relations are understood to be exact in the electron 
and proton masses. The presentation follows [6]. 

Let us introduce the following kinematical A-functions: 

A, ^ A[-(pi + A;)2,-p?,-P] = {yi-y^fS\ 

As = \[-{p^ + h)\-pl-kl] = S^-Am'M^ 

M ^ X{-{p, + k2)\-pl-kl] = {l-yiyS'-Am'M\ 

A, = X[-(p,+Q,)^,-pl-Qf] = yfS'+4M^Ql 

A, = x[-ip,+p,r,-pl-pl] = ylS' + AM'Ql 



where 
and 



1^1 1 — ■97^f^As■; ^1 



\k2\ = OAJVA/' ^2 = 9u(l-^0 



1 A~ /no _ S 



(A.l) 



A. ^ A[-(pi + A)2,-p?,-A2] = {l-y,)'S'-4M\ 

-A^ = r = -{k2 + ky = Z2 + m^ (A.2) 

A(a;, y, z) = x'^ + y'^ + z'^ — 2xy — 2xz — 2yz 
= {x — y — zY — Ayz 
= [x-{^+^zf][x-{^-^zf]. (A.3) 

In the proton rest system, pi = 0, the following relations hold: 

\k\ = oxtVA^, k^ = T^ivi-Vh), 



1 A~ UO — S 



(A.4) 

IP2I = \Qh\ = 2M"^' ^2 = M+^Vh, 

|A| = 2M^^' "^^ = IM'(^~^'*)- 

They may be derived with the aid of the relation [79] 

4M1|pb|' = A[-(pA + PB?, -pI, -p|]|pa=o- (A.5) 

Thus, to any of the momenta in the rest system of the proton corresponds one of the relativistic 
invariant A-functions. From (A.l) and (A.4) it follows that the invariants (3.1) totally fix the 
spatial configuration of momenta of reaction (1.2) which in the proton rest system can be 
drawn as a momentum tetrahedron; see figure 3. 
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A. 2 Kinematics for section 6 

For the calculation of QED corrections to the photoproduction process sets of external and 
integration variables are used which differ from the set yi,Qf,yh,Qh- Further, the infrared 
singularity is treated in a different rest frame. For these reasons, some additional kinematical 
relations are needed here. The presentation follows [6]. 
Let us introduce the following A-functions: 

,2 /^_ I u\2 u2l 



a; = X[-^p,-k2)\-{p2 + k)\-kl] = {S-Qfy-4m'W^ 

y, ^ X[-ip, + hr,-ip, + k)^-kl] = [(l-y,)S + Qf]'-Am^W' 



(A.6) 



2 I /\/f2\2 ^^2ti/2 



{S-W^ + M^Y - 4mW 



With the aid of them and of (A. 5), the following relations may be found in the rest frame of 
the compound system, which is defined by the condition k^ + p§ = 0: 



I^KI ^ _V^ ^ 



2VW^' 2VW' 



2 



' " 2VW^ 2VW^ 



\'^2 I 



^2 ,o,« _ {l-yi)S + Qf 



2VW^' 2VW^ 

The invariant mass of the compound system at rest is W"^ = —{k +^2)^- 

A. 3 Kinematics for section 7 

In the approach of section 7 to the phase space, the first two integrations will be performed 
over the photonic angles in the rest frame R of the (76) compound system: A^ = k^ + k2 = 0. 
For that purpose, the expressions for some four momenta and invariants are needed in that 
system in order to finally express yi,Qf,Zi in (2.40)-(2.42). The velocities which are used in 
the infrared divergent part of the cross section are also calculated in the R frame. 
The necessary A-functions are: 

A- ^ A[-(A-A;)2,-A2,-A;2] = {r-my, 

X, ^ X[-{A-k,f,-A\-kj] = {Ql + z^f + Am'Ql (A.8) 

A2 = X[-{A-k2Y,-A',-kl] = {T-m^, 

and the A^ has been introduced in (A.l). 
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The four-momenta are: 



\" 2 

\tR\ — V fe i-O.R _ T — m _ Z2 



\tR\ _ V^i i,o,i? _ Ql + T + rn^ 

2^ 2v^ (A.9) 

1^1 1 - V^' ^' - 2v/f • 

In the -R system the following cartesian coordinates are chosen: the z-axis be parallel to /cf 
and the yz-plane spanned by the vectors A;f ,pf ; the angle between the 2;-axis and pf be 'dp. 
Then it is: 

k = {A;°'-^sin^i?cosv?R, A;°'-^sin^ijsin v?/?, k^'^cos^R, A;°'^} , (A. 10) 

h = {0, 0, l^fl, A;?'^}, (A.ll) 

^2 = {-A;°'^sin^^cosv?ij, -A;°'-^sin^i^sinv9ij, -A;°'-^cos^ij, A;""^}, (A.12) 

pi = {0, Ipflsin^p, Ipflcos^p, j9?'^}. (A.13) 

From these expressions, one may easily derive (7.9)-(7.12). 

A. 4 Some phase space parameterizations 
A.4.1 The integral d^k2/2k^ 

The following phase space integral has to be calculated in section 3.1: 

d^k2 f \k2 

It is: 

Qf = -2m^ + 2A;°A;° (1 - /^i/^s cos ^12) . (A. 15) 

One has the freedom to identify ^12 = '&■ In the rest system of the proton the formulae of 
appendix A.l may be used. It follows: 

y^ = l-'^kl S = 2Mkl (A.16) 

These two relations allow to replace the energies introduced above by the invariants yi and 
Qj. The corresponding two three-momenta may also be expressed with the aid of A. 4: 

|^J = ^, 1^1 = ^. (A.17) 

' ' 2M' ' ' 2M ^ ' 

Inserting these relations into (A. 14), one gets: 



I ^ = J^-^dk',d COS M^. (A. 14) 



, dyidQl (A.18) 
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A.4.2 The integral dVk 

In section 3.1, tlie following expression has to be simplified: 

dk 



dV, 



2A;0 



5 



{pi + QhY + M^ 6 Ql-{p2~Pi 



(A.19) 



With Qi = Qh + k it is in the proton rest system: 
1 



dVi 



kodko d cos i!}^ dip 6 {pi + Qi-kf + Ml 6 Ql- {Qi-k) 
-kodkodcos^^d<p6[{-M^ + 2pi{Qi - k) + {Qi - kf + Ml] 5 [qI - {Qi - kf 



kodkodcos'd.,dip5[{-M^ + 2p^{Qi -k) + Ql + Ml] S [qI - Q] + 2Qik 
dip dip 



2 

7:kQ 



2 ''2M2\Qi\k^ SM\Qi 



(A.20) 



Here ip is the angle between the planes defined by the three-momenta (/ci, k2) and {k,p2), see 
figure 3. 

Now the integration over dz2 may be introduced with the aid of an auxiliary (5-function 

5{z2 + 2k2ky. 



dT, 



dzo 



mQi 



-d<p6{z2 + 2k2Qi-2k2Qh)- 



(A.21) 



The following coordinate system will be chosen. The z axis is parallel to Qi. The vector /c2 
stays within the xz planeand has an angle i)p with respect to the z axis. The projection of 
the vector Qh = P2 into the xy plane has angle ip to the x axis, and the vector itself has the 
angle -(9 with respect to the z axis. Then, it will be 



k2 ■ Qh = k2 ■ P2 = \k2\\p2\{sm i!}p sin {} cos ip + cosi!}p cos '(9), 



(A.22) 



and 



dT, 



dzo 



lQM\Qi\\k2\\p2\\ sin -(^p sin -(9 sin 99 1 
dz2 



im\Qi ■ ih X P2) 



(A.23) 



The three-momenta used here and in the following may be calculated with the aid of the 

A-functions (A.l). The dTk is related to the volume Vr of the tetrahedron of momenta as 

shown in figure 3: 

1 



Vt = 



6 



Ql ■ ik2 X P2] 



24M 



Rz- 



(A.24) 



Explicitly, 



1/2 



^T= g QfA;2V2-Qf(A:2P2) -A;2^(p2gO -P2(Q^^2)' + 2(QzA;2)(A;2P2)(P2gO • (A.25) 
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The positive-definite function R^ in (A. 24) is tlie Gram determinant of 4- vectors ki, pi, 

kf (kipi) {kik2) {kip2) 



R^ = -Ai{ki,pi,k2,P2) 



-16 



(A.26) 



(A.27) 



(Piki) pI {pik2) {P1P2) 

{k2ki) {k2Pi) kl {k2P2) 

{P2ki) {P2P1) {P2k2) pI 
It may be expressed by tlie G-function introduced in [79]: 

Rz = ~T«T74^(^'' '^'^' '^^' -^9' '^■5'' -^fc)- 

For furtlier use, it is convenient to write Rz in a form wliich exhibits its exphcit dependence 

on Zi or Z2'- 

R, = -A^zl + 2B,zi-Ci = -A2zl + 2B2Z2-C2, (A.28) 

A2 = \ = Ai, (A.29) 

B2 = {2M'QKQ^-Ql) + {l-yi){yiQl-y,Qf)S' + S\l)Qf{m-yH)} 

= - 5i{(l)---(l-y/)}, (A.30) 



C2 = {{l-yi)Ql~Qf[{l)-yH] 



S^ + 4m' 



{yi - VkKviQl - yHQf)S' - m\qI - Q]f} 

(A.31; 



with the discriminant 



Dz = B^2^ ^l,2Cl,2 



1 



64M4 



With (A. 24), the phase space takes the form 

rflfc = 2 X 



dzo 



A(A5, A/, Aq)A(Aq, A/i, Afc) 



dzr>. 



ajk, 2Jr: 



(A.32) 



(A.33) 



where the overall factor of 2 at the right hand side corresponds to 99 G [0,7r] ^h^ z\{2) G 
[2:^2) 5 -^1(2)^] ■ The latter boundaries for the integration over zx(2) are derived in appendix B.l. 

A. 4. 3 The compound phase space volume rfF^g 
The following phase space integral has to be calculated in section 7.1: 

dk2 dk 



dV, 
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2k2 2A;0 



6\A -k2-k) 



dko 



M'^^^-'^^'^- J 2 



ko 



■dkldcosM^5[{K - k2f]. 



(A.34) 



In the rest system of the (76) compound one may use the expressions for the energies and 
momenta (A. 9). In this system, it is 

1 



dkl5[{K - k2 



j dkl5{-T-m^ + 2^k. 



2v^' 



With \k2\ = vA2/2v^ and A2 = A(r, m , 0) = (r — m ) , see appendix A. 3, one gets: 

.2 



r — ra 



ir 



d cos i!}dip. 



(A.35) 



(A.36) 
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AAA The integral d^p2/2p2 



In section 7.1, tlie following integral has to be calculated: 



d%S{pl + Ml)6[{p2-pif-Qt 
2tt I ^^dpld cos M\{p2 - Pif - Q 



6 



iP2 - Pi - hY + r 
iP2-Pi -kiY + T 



In the rest system of the proton it is r = — (p2 ~ Pi — ^i 



Ml 



S 



2Mpl 



(A.37) 



2p2ki. 



In appendix A.l, the expressions for the four-momenta p2 and ki in this frame in terms of 
invariants have been derived. It is — 2p2^i = —{V-^h^s/2M'^) cosi} + . . . This relation allows 
to use the second of the (5-functions for the calculation of the integral over cost?: 



dcos'dd 



-{P2-Pi -kif + T 



2M^ 

2 I /\/f2 I 



With \p2\ = v^/2M and -(p2 - Pi)^ + Ql = -2Mp9, + M'' + M';^ 
allows to perform the integration over dp2 and one gets: 



n 



2^/\^ 



(A.38) 



Ql, the last (5-function 



(A.39) 



A. 4. 5 Notations 

Here, we confront some old and new notations. This may prove useful if someone tries to read 
the preprints [4], [5]: 



S 



X 



Si -- 


= yiS, 


T = 


= VhS, 


t = 


= Ql 


Q' - 


= Ql 


X = 


= {i-yi)s, 


Y = 


- Ql 


S-X -- 


= yiS. 











(A.40) 



B Kinematic boundaries 

We now come to the derivation of the many different kinematical boundaries which are con- 
nected with the different choices of the integration variables and their sequential order. By 



inax(min) 



f"'"'"^"""'(Vi, V2, . . .), the maximum (minimum) value off will be denoted which may be taken 
at fixed values of Vi, V2, . . . Further, the extreme value of variable v is v\ it depends only on S 
and possibly on M, m. 

The physical regions of invariants (3.1) may be found from the conditions of existence of 
the momentum tetrahedron depicted in figure 3. This tetrahedron is defined by the triangles 
A(A;i, — A;2, —Qi), ^{Qu ^P2, —k), and A(A;i, —A, — P2) in the proton rest system. The areas of 
these triangles are equal to 

y-xiQipv^p), 



y-\{kikiQf), 

correspondingly; they exist if 



-xikiA^pr 



KklklQl < 0, \{Qlp2ie) < 0, \{klKlp2^) < 0. 



(B.i; 
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With the expressions (A.l) for the squared three-momenta in terms of A-functions, the follow- 
ing boundary conditions are derived: 

A(A5,A,,A,)<0, X{Xg,Xh,Xk)<0, A(As, A,, A^) < 0. (B.2) 

The boundary equation for the first condition is: 

M^Qt + QfyiS^ + m'yfS' - XsQf = 0. (B.3) 

The boundary equation for the second condition is: 

yfS'Ql + ylS'Qf - M\Q^ - Qlf - yMQl + Ql)S' = 0. (B.4) 

It is symmetric with respect to the interchanges {yi ^^ yh) and {Qf ^^ Ql) . 
The third boundary condition may be rewritten: 

i2M^Z2 + yhS^ + 2M^Qiy - XsXh < 0. (B.5) 



The Z2 is known to be positive definite (see appendix B.1.1). The maximal value z1^^^{yh, Q 
of Z2 at given values of {yh, Ql) is one of the roots zf{yh, Ql) of (B.5), 

1 



4{yh,Ql) 



iyhS' + 2M'Ql)±JXsX 



(B.6) 



2M2 
namely ^2'. It is positive definite if 

M'Qt + QlyuS' + mVhS' - XsQl < 0. (B.7) 

From the zero of (B.7) the third boundary equation results. By the replacements Q\ <-^ Q^ 
and yi ^^ yh it transforms into the boundary equation of the first condition. 

B.l Boundaries for zi, Z2 

B.1.1 Boundaries for zii2){yuQf,yh,Ql) 

Here we closely follow [6]. From (A. 32) and (B.2) one may conclude Dz > 0. The invariants 
zi and Z2 are positive definite. This may be understood easiest from the general statement 
that an invariant of the form I = —2pk is positive definite when the momenta p, k are on the 
mass shell. One may also use an argument based on the photon energy in (A. 9), which is 
proportional to ^2. For Zi the same property may be derived by a replacement of k2 in favour 
of ki in the arguments which lead to the lower boundary of Z2 (i.e. in section A. 3). 

From (A.l) and (A. 29) it follows Ai^2 > 0. The Rz which is introduced in (A. 24) and 
rewritten in (A. 28) is proportional to the volume of the tetrahedron in figure 3, i.e. it is also 
positive definite. This volumee vanishes at the boundary values of the physical regions of the 
invariants. Rewriting Rz, 



Rz = Xq 



^1(2) ■" ^1(2) ^1(2) - ^1(2) 



(B.8) 



one may derive an expression for the boundary values of zi and Z2 at fixed values of yi, Qf, 
Vh, Ql- 

max,min/ ^2. n'i\ — -^1.2 =1= V -Pz _ W,2 



max, mm/ /^2 /02\ J^i^ v -" z ^i^ /-n r\\ 

%,2 {.Vh Qf; yh, Qi) = -. = - — —y^- (B.9) 

Az -Dl,2 =F V-LJz 
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The Dz in the above relation must be positive. Since the Zi(2) are positive definite, the -Bi(2) 
must be positive, too. Finally, from the second equality in the above chain, the Ci(2) are also 
seen to be positive. In the ultra-relativistic approximation, the latter is also evident from 
(A.31). 

We mention here that the integration limits for Z2 at fixed values of yh-, Q\ will be derived 
in appendix B.4 and at fixed yjB, Qjb in appendix B.5. They are also infiuenced by the 
photonic cuts, which are introduced for leptonic variables in section 6; see appendix B.1.2. 

B.1.2 Boundaries for zi(2){yuQ'i,yh,Q'h) with photonic cuts 

In section 5 we argued that cuts may be applied on the momentum of the bremsstrahlung 
photon. This is of some interest in the case of leptonic variables in the region of small x and 
large y. In this appendix, the necessary formulae will be derived. 
We define the photonic cuts by the following conditions: 

COs(^--)=C^i„>COS^^ > C^ax = COs(^--^). ^ ' "^ ^ 

With these cuts, photons with an energy larger than R ™"i and within a cone defined by 
cos(^ ™^^'™°) are excluded from the observed cross section. 

Both experiments at HERA have a 7 counter which may register photons with an energy 
E^>1 GeV and an emission angle < ^^ < 0.5 mrad in the HERA laboratory system. 

In order to incorporate the cuts (B.IO) into the cross section calculations one needs two 
invariants which depend exclusively on E^ and 9^ and, possibly, on Ef. and Ep^ the energies of 
the electron and proton beams in the laboratory system. 

One candidate is zi. 

zi = -2kik = 2EeE^{l-f3eCose^), (B.ll) 



m^ 



Another one may be found to be 

I^ = -2p,k = {yi-yh)S = 2EpE^(l + f3pCose^), (B.13) 



f3p 



M2 
.1-^. (B,14) 



Dividing (B.ll) by (B.13) one gets a first condition which is formulated completely in terms 
of integration variables and measurable parameters. A second one may be obtained by elimi- 
nating the cos^^ from (B.ll) with the aid of (B.13): 

mill _ ^e[^ ~ PeCminj c( \ ^ ~ ^ -^ejl ~ PeCmaxj c( _ „, \ — ^max 

^c = p (T _i_ fi r W^i'' ~yh) ^ Zi ^ P (T \ fi r W^i'' Vh) — Z^ , 
P\ ' h'p^vam) P\ ' r'p'-^nisx) 



zf- = 2E,EUl + ^)-^^S{yi-yh) < z,. 



(B.15) 
Taking into account Z2 = Zi + Qf — Q\ one may realize that (B.15) represents cut conditions 
on the first integration variable. 
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The actual modifications of tlie analytical integrations over zi(^2) due to the cuts may be 
found in D.l. After this first integration, the finite hard cross section part undergoes two nu- 
merical integrations, which are not essentially affected by the cuts. Concerning the treatment 
of the infrared singularity, nothing changes compared to the presentation in section 4. The 
reason is simple. We subtracted from and added to the complete bremsstrahlung correction a 
simplified expression, well adapted to the infrared problem. Now, in the presence of cuts, we 
have the freedom to subtract and add the same expressions. Here, it is of importance that the 
applied cuts do not influence the infrared behaviour since the contribution of hard photons is 
manipulated. 

Concluding, 

zr = max {zrivi, Qf, Vh, Ql), zf^{E,, E,, c^i„, yi, Vh), ^^ (^e, E^, m, y^)} , 



B.1.3 Boundaries for ^2(1/^, Q^) 

The following boundaries are used for the calculation of cross sections in hadronic and Jaquet- 
Blondel variables where the integration over r is the third one. After two integrations over 
photonic angles the ultra- relativistic approximations are justified and will be applied. For 
given values of yh and Q\, the upper limit of Z2 is given in (B.6) while the lower one is zero; 
compare the expression for k^'^ in (A. 9): 

0<Z2<Xhil-yh)S. (B.17) 

The boundary values of Z2 define also those for the variable r when the latter is used as the 
argument of the last integration at given values of either yh, Q\ or yjB, Qjb- 

B.2 Boundaries for leptonic variables 

From now on, we will use the following notations for external variables: if a variable is maximal 
or minimal as a function of the other one (and of S and the masses), the superscripts min or 
max will be used to indicate this. Absolute extremes which depend only on S and the masses 
get a bar. 

The maxima and minima of an internal variable at fixed values of the other internal variable 
and of the external variables get the subscripts I, III, and II, IV, respectively. In the figures 
these functions are boundaries. If the integration order is changed the corresponding curves 
in the figures get a different analytical meaning. This will become clear in the following 
discussions of integration regions. 

B.2.1 External variables £1 = {yi,Qf) 

We start with the derivation of absolute bounds for W^. 
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For the invariant hadronic mass M| the lower hmit is defined by the kinematical threshold 
of hadron production, i.e. the mass of the lightest strongly interacting particle, the pion. The 
upper limit of M|, for a given W^, may be derived from the expression for /cq in the rest 
system (A. 7). It is reached for k^ = 0: 

{M + m^f <Ml<W^. (B.18) 

The above inequality may also be interpreted as the absolute lower limit of W'^. 

The upper bound for W"^ = — (pi + ki — k2Y follows from a consideration of the defining 
four momentum squared in the centre-of-mass system: 

(M + m^f <W'^ <W'^ = {^rs- mf = {VS + M^ + m^ - mf . (B.19) 

One may eliminate W'^ from (B.19) with the aid of (6.2)®: 

2 



^<yiS-Q'i < NS + M2 + m2 - m) -Ml (B.20) 



The boundary equation (B.3) originally depends on the variables we are interested in: yi and 
Q;^. It may be pursued once more. The boundaries for yi at a given value of Q^ shall be 
determined. The lower bound of yi results from (B.20); the upper limit is the solution of 
(B.3): 

yrm = ^<yi< yr\Ql). (B.21) 

where 



yr^iQl) 



2171^ 



-^y^s^m — Q 



< m. (B.22) 



Here the yi is the maximum of yi and may be determined by differentiation of (B.22) with 
respect to Qf: 

S - 2mM 
m = — ^ — • (B.23) 

The absolute minimum of Q^ is zero. The condition |/™"(Q^) = vT^^iQf) is fulfilled at 
Q^ = and also at the absolute maximum of Q^. It may be used for a determination of the 
limits of the latter variable: 

The physical region Si = {Qf,yi) is shown in figure 33. 

B.2.2 External variables S'l = {W^,Qf) 

The bounds for W'^, (B.19), have been derived in the foregoing appendix. 

Those for Qf at a given value of W'^ may be obtained from (B.3) by replacing the yi by 
W^ with the aid of (6.2): 



From the lower limit in (B.25) the estimation of the lowest accessible Q^ (6.3) in the photo- 
production region was derived. 

The physical region £'i = {W'^, Qf) is shown in figure 34. 



^If not stated differently, the pion mass is neglected here and in the following. 
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Figure 33: Physical region (Qf, yi) for the cross section in leptonic variables. 




(M + mj' 



Vi^ W^ 



Figure 34: Physical region (W , Qj) for the photoproduction process. 
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B.2.3 Integration variables Xi = {yh,Q 



For the hard photon corrections in leptonic variables, the last two integrations are numerically 
performed over Xi = {i/h, Ql)- The hard part of the infrared divergent contribution is integrated 
analytically; thereby the last integration is that over yh- 

Inserting (6.5) into (B.18), one gets a first condition the left hand side of which is equivalent 
to Xh = 1: 

0<yhS-Ql<W^~M\ (B.26) 

The (B.4) defines a pair of intersecting straight lines : 

Qkn = Qf + ^{yi - yk){yiS ± ^/\). (B.27) 

From the second relation, the lower bound of Ql is obtained: 

Qr''{yk,yi,Q^) = Qlu- (B.28) 

The upper bound of Qj^ depends on the actual value of yh- The two curves Q\ = yhS and 
Q\ = Q\^ meet at the maximum value of Q^ which is reached at the following value of yh'- 

r = M^S- (B.30) 

At yh > yhd it is Qf^'^^'iyh, yi, Qf) = Qli- 

The minimum of yh may be obtained from the equality Qln = yhS; see figure 35: 



yi(yiS-J\) + 2Qfr 
yr\yu Ql) = -^ ^ ^- (B.31) 



The maximum of yh at a fixed value of yi is equal to yi and independent of Qf; this may 
be seen from the expression for k^ in (A. 4). The infrared singularity is located at J-'{yi,Qf)- 
When calculating the hard part (4.13) of the infrared divergent cross section contribution, one 
has to take into account an additional condition: 

A;° > e > 0. (B.32) 

In the approach to the phase space where the boundaries of this appendix are used the infrared 
problem is treated in the proton rest frame. From the expression for k^ in (A. 4) the following 
modification of the integration boundary results: 

y'r = yi^ y'n^) = ^' - ^- (B-33) 

The net physical region for X/ = {yh, Ql) is: 

Qlii < Ql < yhS if yh < yh^, 
Qlu < Ql < Qli if yh > yn,, (b.34) 

y^iyuQl) < yh < yr\^- 
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Q^ 








F(y^Q;') 



yh Yhi 



YhZ Jhd 



Jh 



Figure 35: Integration region (y^, Q\) for the cross section in leptonic variables. 

When calculating the twofold integrals in (4.13) over the kinematical domain, which is 
derived here, one will need the substitutions of \/Ci and \/C2 at the kinematical boundaries. 
The square roots at (QDi.n '^^^ t)e taken exactly: 



CiiQ 



C2{Q' 



hi) 



C,iQ' 



hll) 



C2{Q 



hll) 



{yi - yh)S 
{yi - yh)S 
{yi-yh)S 
{yi-yh)S 



S 



Ag + yiS 
2M2 



+ Qf 



K + yiS 



^ 2M2 ^' 



K — yiS 



(B.35) 
(B.36) 
(B.37) 
(B.38) 



The corresponding roots at Q^ = ynS may be taken in the ultra-relativistic approximation in 
the electron mass; see the expressions (A. 31). One arrives at the following absolute values: 

(B.39) 
(B.40) 



'Ci ~ S'^\yh-xiyi{l+yh-yi)\, 
^C2 - S^\yh{l - yi) - xiyi{l - yh)\, 
which vanish inside the region of integration over y^, (B.34), at 

xiyi{l - yi) 



yhi 



1 - yixi 



[B.41) 
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correspondingly. Within the kinematical hmits of figure 35, the following inequalities hold: 

yr < Vh, < yH2 < VK, < y'n^). (B.43) 

The integration region Xi = {y^, Ql) is shown in figure 35. 

B.2.4 Integration variables J^' = {Ml,QD 

The boundaries for M| have been derived in (B.18): 

{M + m^f <Ml<W^. 

The infrared singularity is located at J-'{W'^,Qf). When calculating the hard part (4.13) of 
the infrared divergent cross section contribution, one has to take into account the additional 
condition for the photon energy (A. 7): 

^o,« ^ W -Ml ^^^^ 

This condition may be transformed into the following modification of the upper bound of M|: 

(M + m^f <Ml<W^- 2eVw^. (B.45) 

The limits of the variable Q^ for given values of W'^,Ml,Qf will be derived now. For that 
purpose, one has to replace in (B.4) the yi with (6.2) by W"^ and Q^ and the yh with (6.5) by 
Ml and Q\ and solve the resulting inequality for Q\ in terms of Ml, Qf, W"^ [7]: 

Qhu = g^^^^■^^"V^ Ml Qf) = ^ [{W - Ml) {mS ±^,)+ 2MlQf] . (B.46) 

The two straight lines meet in the infrared point while the two extreme values of Q\ at given 
external variables are reached for Ml = (M + 171.,^)'^. 

The physical region 2^ = (M|, Ql) is shown in figure 36. 



B.2.5 Integration variables X^ = {xh,Qi) 

For the discussion of the Compton peak we need the integration boundaries for the variables 

n 



Xh and Ql^ where the latter has to be integrated over first. 



Using in (B.4) the relation y^ = Ql/{xhS), one immediately arrives at the boundaries for 
Ql^ at given values of Xh, yi, Qf'- 



Qkn = Ql \yuQlx,) = 2[x,y,S-Q^,l(x,MY] ' ^""-''^ 

At the lower limit of Xh the two curves Qfim meet; see figure 37. From (B.47) one derives 
that this is realized if Xh = xi, which is at the same time the infrared point. The upper limit 
for Xh may be got by rewriting (6.5) and inserting (B.18): 

Ql/^h -Ql = Ml -M^> m^{2M + m^). (B.48) 




(M + rnj' 



M;,^ 



Figure 36: Integration region (M^,Ql) for the cross section in leptonic variables. 



Q. 



F(X(,Q,') 








1 X/, 



Figure 37: Integration region (xh, Q\) for the cross section in leptonic variables. 
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Neglecting here the pion mass, the absolute maximum of Xh is evidently equal to one: 

xi<Xh< 1, (B.49) 

Taking the pion mass into account, the value Xh = I is forbidden; in fact Xh = I corresponds 
to the kinematics of the elastic channel. 

The integration region Xi = {yh, Ql) is shown in figure 37. 



B.3 Boundaries for mixed variables 

B.3.1 External variables Sm = {VhiQf) 

The boundaries for Qf are the same as in the case of leptonic variables; see (B.24): 



< g/ < 



5 + M2 + m2 



= Q^ 



The maximum of yh in the two-dimensional distribution [yh{Qj)]yi{yh,Q'h,Qj)\ is, inde- 
pendent of the value of Q^, equal to yi. Thus, its absolute maximum is given by the value for 
yT^^iQf) of (B.22). The minimal value of yhiQf) niay be obtained by studying the functional 
dependence of y^h^^^iyuQf) in (B.34) on yi. The extremum is reached at yi = yf^'^iQf)- The 
explicit expressions are: 



yr\Qh < Vh < yr^'iQf), 



(B.50) 



where 



yrm 



QUJiyr^'Sy + ^M^Qt-yTS 



{jjf^^-Sy + 4M2g2 + ^max^ _ ggS 



yr^'iQl 



s 



1 



g? (v^ + Qf) {s - v^) 



< 



01 

S'' 



2m2 



— yA^Am — g 



< 



yi- 



At the right hand side of (B.51), the yf'''^ stands for yf'^'^iQ'j) [see (B.22)] and 

A™ = gf (g? + 4m2) . 

The second part of (B.51) has been obtained from the identity 

1 



[B.51) 
[B.52) 

[B.53) 



(^max^)2 + 4JVf2g 



2^2 > *^V ^- 



QfVA 



The two curves yf^^^iQf) and y^"^{Qf) meet each other at Qf = Qf as defined in (B.24). 
With the variables Sm one may define a corresponding set of scaling variables: 



*Vm ^l ? 



yh, 






-2piQh yhS' 



(B.54) 
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Figure 38: Physical region (yh, Qf) for the cross section in mixed variables; lower part for 
Xm > 1, upper part for x^ < 1. 

Here the mixed variable Xm is defined from botli leptonic and liadronic variables, while 
ym is defined exclusively by the hadronic kinematics. For y^S > Qf, the scaling variable Xm 
fulfills < Xm < 1, while at yhS < Qf it is Xm > i- 

Values of x which are larger than 1 seem to be exceptional. Nevertheless they are physical 
in the case of Xm- This statement may be checked e.g. by reproducing the phase space 
volume (C.7) in mixed variables. 

The physical region Sm = {Qf,yh) is shown in figure 38. 

In the kinematical region with Xm > 1 it is ensured that Qf — Q\ > Q since yhS = 
Q\ + M^ — M^ < Qf. Thus, the infrared point may not be reached and only the radiative 
process is possible. The equality in (B.51) is realized for yf^^S = Qf; then it is y^^'^S = Qf. 

B.3.2 Integration variables X^ = (VhQh) 

The complete analytical integration of the hard part of the infrared divergent cross section 
will be performed first over Qf and then over yi. We consider only the case Xm < 1 explicitly. 
As in the case of leptonic variables, the bounds for Qf at given values of yi, yh, Qf niay be 
derived from (B.26) and (B.27), and figure 39. The Qf^^^iyuVhiQf) is either equal to Qf^ or 
to yhS; both curves meet at the value yi^: 



yid 



1 



l+Xm + {l 




[B.55) 



with r = M'^/S. 



84 



The lower bound Q'i^"\yi,yh,Q'i) is equal to Qln- 

S I — 

Qhll = Qf + j^iyi-yh){yiS±^\). 

The upper bound for yi is dependent only on Qf; see (B.22): 

1 



y-riQl) = 5^ 



y ^s\n ~ Ql 



The lower bound for yi is reached at the infrared singularity J-'{yh, Qf)- In the calculation of 
the hard part (4.13) of the infrared divergent cross section contribution one has to take into 
account an additional condition: 



r > e > 0. 



[B.56) 



The infrared problem is treated in the proton rest frame. Corresponding to (A. 4), this condi- 
tion may be transformed into the following modification of the integration boundaries (B.34): 



Vi 



yh - yr\e) =y^ + 



2eM 



(B.57) 



The integration boundaries are: 



Qlu 


< 


Ql 


< 


yhS 


if yi > yid: 


Qlii 


< 


Ql 


< 


Qli 


if yi < yid, 


r(e) 


< 


yi 


< 


yr'-m- 





fB.58) 



The integration region Xm = {yi, Ql) is shown in figures 39 and 40. 

In appendix D.3.3, the square roots of the functions Ci and C2 will be needed. In ap- 
pendix B.2.3, we discussed these roots at the integration boundaries for the case of leptonic 
variables. Here, exactly the same discussion applies with inclusion of the relations (B.35)- 
(B.40). Although, there is one modification which concerns the location of the zeros of Ci and 
C2 at the boundary Ql = ynS. Differing from the leptonic case here they are at some values 
of yi for fixed external values of yh and Qf: 



yii = i + yh , 



yi2 



l-a;„(l -yh)- 



They satisfy the inequalities 



yii <yh< yid < yi2 < yT 



(B.59) 



(B.60) 



if additionally the condition yh < 1/(1 + r) is fulfilled. In this case, the yn does not belong to 
the physical interval {yh, yi"^""^) of variation of yi, see figure 39. 

We only mention that one needn't know the values of the roots of Ci and C2 on the fourth 
boundary since the integration is performed along its direction in the phase space. 
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yu yid yz2 yi Ji 



Figure 39: Integration region (yi,Q\) for the cross section in mixed variables, Xm < 1- 
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F(y^,Qi^ 







yh^ 



Q hll 
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Figure 40: Integration region (yi, Q\) for the cross section in mixed variables, x^ > 1. 
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B.3.3 Integration variables 2^ = {Q\,yi) 

Here the boundaries for the integration of the hard, infrared finite part of the bremsstrahlung 
cross section are derived. In this case, the integrals over yi are performed analytically in 
appendix D.4. The last one over Q\ has to be done numerically. 

The lower limit of yi at fixed values of Qf, yh, and Ql may be found from (B.4): 

(B.61) 



2Q 



yrivk, Ql Qf) = 7^ VkiQl + Ql) + -^\Ql - QlW^ 



s' 



There is only one solution if Xm > 1 since then it is definitely Qf > Ql- For the more common 
case Xm < 1 the solution has two branches. Geometrically, they coincide with the curves 
Qhiii in figure 39. The two curves (B.61) meet at the infrared point J-'{yi,QD. This point 
marks the lowest allowed value for yi since yi > y^, see (A. 4). There is no need to exclude 
the infrared singularity from the integration boundary, since the integrand has been regulated 
there. The upper limit yf^^iyh, Qh Qf) of Hi does not depend on yh or Q\ and is as derived 
earlier [see (B.22)]: 



yr^'iQ'i 



2rn? 



^S^m ~ Ql 



<yi- 



yr^'S-Jiyr^SY + ^M^Q'i 



The upper limit of Q\ derives from the inequality M| > M^ together with (6.5). It is 
independent of the actual value of Qf: Ql^^'^i^yh, Qf) = yhS. The curves y^^^{yh, Ql,, Qf) and 
yr^^iQl) iiieet at Ql = <5|™°; this condition may be used for the calculation of (5|™"^(|//i, Q^)'- 

Qr^y^m = Ql + i^Avr"" - yn) 

S {yr"" - yh) (Qf + v^) 
= Qf- -—^ ^^li \ (B.62) 

Taking all the above conditions together, the physical region for the integration variables 
is found: 

Qr\yk,Qf) < Ql < yhS, 

(B.63) 

yriyH,QlQl) < yi < yr''iyh,QlQl)- 

For the analytical integration over yi in appendix D.4, we need the expressions for the 
roots of Ci and C2 at the integration boundaries. 
At y™'^'', they are: 

^;(^ ^ ^s[Qt + 2m\Q} + QX)]- yX^SjQl + 2m'yh) ^ ^^^^^ 

^^(^ ^ V^sjQlQl + 2m\Qf + Ql)] - V)^S{Ql + 2m^y,) ^ ^^^^^ 

Since the lower boundary of yi has two branches, the roots have different analytical ex- 
pressions there and the integration region has to be divided into two subregions: 
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• Region I: Qf > Ql 






. Region II: Qf < Ql 






The following relations hold: 






^c^ivr)' = 


{Qf-Ql)[2SQl-QKyhS + VyH)] 


(B.67) 


2QI 


^C2{yf') = 


{Qf-Ql)[2SQl-Ql{y,S-VK)] 


(B.68) 


2QI 


^\{yr)' = 


ykS{Qf-Ql) + {Q^ + Ql)VK 


(B.69) 


v/ci(yr)" = 


iQl-Qf) 2SQl-QKyhS-v%:) 


(B.70) 


2QI 


v/ci(yr)" = 


(Ql-Qf)[2SQl-Ql{y,S + y%:)] 


CR 71 "l 


2QI 


[D.l i) 


A(?/r)" = 


y.SiQl-QD + iQ^ + QDVXi: 

o/n2 


(B.72) 



All these roots are positive definite. Again, the values of the roots at the fourth boundary are 
not needed since the integration is performed along its direction. 



B.4 Boundaries for hadronic variables 

B.4.1 External variables Sh = {yh,Qh) 

The physical region of Sh = {yh, Ql) ^^^J be derived from the boundary condition (B.7) and 
from the inequality 

0<yhS-Ql = W^-M^ <W^- M^ (B.73) 

where W^ may be taken from (B.19). This proceeds exactly as in the case of the leptonic 
variables, by a replacement of {yuQf) by {yh,Qh) everywhere in the argumentations. 
One gets for £h the following physical region: 



< Ql < 






m 



Ql 



yr\Ql) = ^ < yh < yr\Ql). 



where 



yr^'iQl) 



2m? 



-^\J\sQl{Ql + 2m^)-Ql 



<yh- 



Here, the yh is the maximal possible value of yh'- 

S - 2mM 

yh = ■ 



s 



(B.74) 



[B.75) 



(B.76) 



The physical region Sh = {Ql,yh) coincides with that for the case of leptonic variables, 
which is shown in figure 33. 



B.4.2 Integration variables Xy, = {yi,Qf) 

Similar to the case of the mixed variables, for the determination of the hard part of the infrared 
divergent contribution we first integrate over Qf and afterwards over yi. This will be done in 
appendix D.3.4. 

Two bounds for Qf at given values oiyi, yh, Ql may be derived from (B.27) by interchanging 
yi with yh and Qf with Q^ as was indicated in connection with (B.4), 



Qi i,ii 

and two others from (B.3), 



Ql + 



s 



2M2 



{yi - yh) l±y>^h - yhS) 



Qf 



IIIJV 



2M2 



'1 



m)S^-Am'M'±,/x^i 



(B.77) 



fB.78) 



As may be seen in figure 41, the absolute maximum yi of yi divides the curved boundary into 
two parts. The yi is defined in (B.23). It is the root of A^: \i{yi) = 0. Thus, the two curves 
Qfiii IV iiieet there. Further, the minimum of yi at given y^ is equal to yh and independent of 

When calculating the hard part (4.13) of the infrared divergent cross section contribution, 
one has to take into account an additional condition: 



r > e > 0. 



[B.79) 



The infrared problem is treated in the proton rest frame. With the expression for k^ in (A. 4), 
this condition may be transformed into the following modification of the lower bound for yi: 



yl 



yh 



yr\e) 



yh 



2eM 



(B.80) 



For a complete definition of the integration region one needs yet the points where each of the 
two straight boundaries meets one of the two curved. The points yi^ and 1//2 are determined 
from the conditions Qf^ = Qfjy and Qfj = Qfm, respectively: 

[(1 ~yh)iS±V>^ + 2yhS)S - 2M\Ql + 2m')]{^h T yhS) 



yi2{3) 



T 



2S[{1 - yh^Wh T yhS)S ± 2M^{Ql + m^)] 
2M'Ql[il-2yh)STV>^ 



2S[{1 - t/^)(v^ T yhS)S ± 2M^{Ql + m^)] 
In the ultra-relativistic approximation. 



yi2 



yi3 



2S + yh{^h-yhS)-2rQl 
2S + v^ - yhS 



yh 



^rQl 



yhS + vXh 

The integration limits may be summarized as follows: 



Qf'^'\yi,yh,Ql) - 


- Qh 


if yr°(e) < 


yi 


<yi3, 


QT'^'iyuyh^Ql) = 


= Qiw 


if yi3 < 


yi 


<yh 


Qf^^'^'^hjuyh^Ql) -- 


= Qh 


if yr{e) < 


yi 


<yi2: 


Qf^^^^^iyhyh^Ql) = 


= Qi III 


if yi2 < 


yi 


<yi. 



(B.81; 



fB.82) 



[B.83) 



[B.84) 
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Figure 41: Integration region (yi, Qf) for the cross section in hadronic variables. 



The integration region Xh = {yijQf) is shown in figure 41. The infrared divergence is 
located at J-'{yh,Q'i)- 

In the course of the integrations, the values of the roots of Ci and C2 are needed at the 
boundaries. They may be determined exactly: 



^liQ 11,11) 



C2{Q'lui) 



CiiQf 



111,1V ^ 



C2{Qf 



UIJV) 



{yi - yh)S 

{yi - yh)S 
1 



[l-yi + yh)S- 






[i-yh)S- 



2M2 
^iTyhS 2 



2M2 



2M2 



'l-yi + yH){l - yi)S' - 2M\Ql + 2m 



1 



2 S/fS 



'l-yi + yH)S'-Am'M 



2M2 



2/\/r2 



[l-yi){l-yh)S'-Am'M 



±J\i {l-yh)S'-2M\Ql + 2m 



(B.85) 
(B.86) 



[B.87) 



(B.88) 



B.4.3 Integration variables 2^ = (Qf, 



yi] 



Here, the boundaries for the integration over leptonic variables will be derived for the case 
where the integration over yi is performed before that over Qf. This is needed in appendix D.4 
and will be used for the calculation of the finite hard part of the cross section. We remind that 
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in this case the boundaries for yi are the same as in the case of the mixed variables. Thus, 
the minimal value of yi for given values of Qf, yh, and Qf^ is given by (B.61), 



yr{yk,QlQl 



2^ - -L_ 
2 



2Q. 
and the maximal value by (B.22): 

1 



VkiQf + Ql) + ^\Q1 - Qll^/^H 



y'riyk,QlQf)- ^ 



2m? 



C, V ^S^m - Qf 



<yi- 



For fixed y^, Ql, one may see in figure 41 that yf^^^ and y™™ are equal for the extreme values 
of Qf. Solving the resulting equation, one gets: 



w ^ Qt[STV\s][STV\s-yhSTV\h, 

2S(y,STV^){Ql + m^y,)+4M^QUQl + m') 



Concluding, the physical region of X^ = {Qf,yi) is: 

Qr\yH,Ql) < Qf < Qr\yH,Ql), 
yriQhyh,Ql) < yi < yr^iQhyH,Ql)- 

B.4.4 Variables in section 7: Sh = {yh^Qfi)-, ^h = (t) 

In this section, we restrict ourselves to the ultra-relativistic approximation. 
The physical region for the external variables is as derived in section B.4.1: 

< y?. < 1, 

Q < Ql < yhS. 



(B.90) 



(B.91) 



The limits for r = Z2 + rr? at given values of yh and Q\ follow from those for Z2 which 
are derived in appendix B.1.3. For the case of the hard photon contribution to the infrared 
singular cross section part one has also to take into account the condition 

A;° > e > 0. (B.92) 

The infrared problem is treated in the rest frame of the (76) compound system. There, this 
condition may be fulfilled by the demand 

Z2>Z2 = 2me. (B.93) 

The boundaries are then: 

m^ + 2me<r<m^ + Xh{l - yh)S ^ x^il - yh)S. (B.94) 

At the lower boundary, one has to retain the infiuence of the electron mass in order to ensure 
the positive definiteness of the phase space measure and to prevent unphysical singularities 
during the integrations. 
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B.5 Boundaries for Jaquet-Blondel variables: Sjb = (z/jb^Qjb)? 



In this section, we restrict ourselves to the ultra-relativistic approximation. The definitions of 
the Jaquet-Blondel variables are: 

ym = Vh, (B.95) 

Q'jB = ^{V2L?- (B.96) 

The normalization M has to be adjusted such that in the limiting case of the Born kinematics 
it will hold QjB = Q\- We will immediately see that 

M=- ^^, (B.97) 

l-yh + ■■■ 

where the dots stand for some terms proportional to m^ and M^. An exact expression for 
QjB in terms of yh-, QIj t iiiay be obtained in the rest frame of the proton. In the momentum 
tetrahedron of figure 3 it is /ci — p2 = A. In this system, these three three-vectors span a 
triangle whose area may be expressed in two ways: 



l\ki\\p2±\ = ^-\\{h,P2\A'). (B.98) 

The three three- vectors squared are expressed by invariants in (A. 4): 

iP2±f = -^]^^(^^' >^h, K). (B.99) 

With this expression, the Qjb is determined after the normalization has been adjusted such 
that for Z2 = the Qjb becomes Qj^: 

If the electron mass is neglected the Qjg becomes 

r,2 (1 - yh)Ql - yhZ2 - {Ql + Z2fr/S mmn 

"^^^ = ii-y.)-Qlr/s ' ^""-'''^ 

where it is r = M'^/S. Neglecting also the proton mass, one finally obtains 

Q% = Ql--r^Z2. (B.102) 

The hadronic and Jaquet-Blondel-transferred momentum squares are related by yh = yjB 
and the additional variable Z2. This property may be used for a derivation of the kine- 
matical boundaries. We start from the hadronic integration limits which were derived in 
appendix B.4.4. A chain of simple changes of the integration region allows a derivation of the 



92 



boundary values for Jaquet-Blondel variables: 

3 1 (i-j/h)Q|/j/h 1 VhS {'^'yh)Ql/yh 

Ts = / dQl J dyh / dz2 = dyn / dQl / dz2 

Ql/s 

1 (l-J/fc)-? j/hS 1 {l-yh)S yhS-yhZ2/(l-yh) 

dyh / dz2 / (iQl = / rfyjB / (i22 / C^QjB 



yhZ2/ii~yh) 

1 yhS {l-a;jB){l-J/jB)S 5 1 {l-a;jB)(l-J/jB)S' 

/ c?i/jB / dQj-Q / dz2 = dQl^ j dy^B j dz2. 

Qljs 



(B.103) 



Thus, we derived: 



12 



QjB 



0< g^ <5, (B.104) 

< t/jB <1, (B.105) 



S 

m^ + 2mt< T < m^ + (1 - xjb)(1 - t/jB)*? 

^(l-xjB)(l-yjB)5. (B.106) 

Here we took into account that the lower bound for the r integration is modified in case of the 
hard contribution to the infrared singular part; see the corresponding remarks in the foregoing 
appendix. 

The physical region £^jb = (^jb^^jb) is shown in figure 42. 

C The phase space volume 

In the above appendices, many boundaries for external variable sets and for integration vari- 
ables have been introduced in different combinations of variables and in different order. A 
useful tool for numerical checks of the correctness of all these boundaries is the integral over 
the phase space volume, which must be identical for all the different presentations. 

The phase space volume (3.9) is calculated now in terms of the variables which were 
introduced for the study of the photoproduction process: 

dyidQfdyhdQl- ^^^^ 



4^XsJ ' WRz 

^//T7,mA2 „/n2 max „TJ/2 „/n2 max „ -.max ; 



dW' l^_ dQf _ dMtl_dQl_^^. (C.l) 



A^/Xs J{M+m^)'^ jQ'f'^'"' J{M+m^Y JqI"'''^ J z^^"^ ^/R 
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Figure 42: Physical region (Q'^^jVjb) for the cross section in Jaquet-Blondel variables. 



The boundaries for Z2 are (B.9), those for Q\ are (B.46), and those for Q'f are (B.25). We first 
apply (D.2). Next, the integration over Q\ yields 



dQl = Qr-" - Q 



h ? 



(C.2) 



which together with the result of the first integration and with (A. 29), A2 = Xq, leads to 
fortunate cancellations: 



n 






dW'^ dQf dMl- 



The two subsequent integrals over M| and Q^ are also trivial; after neglecting m.,^: 

t2 i\j2\2 



r 



si 



The last integration yields 
8S\ 6 



dW^^ ' ^{s -W^ + M2)2 - Am^W^. 



5^ + lOm^s - 5M^s + m* - 5m^M^ - 2M^ 



1 



m'^s{s + m^ - 2M^) + -(s + m'^)M^ 



In 



s + m^ - M^ + ^/X 



+ -(s-m2)MMn 



S + 171"^ - M^ - ^/X^ 

(^s-m^ + M^ + v^A^) (^s-m^ - M^ + ^/X^^ 
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(C.3) 



(C.4) 



(C.5) 



with s = S + rn? + M^. Neglecting the electron mass m one gets 



r ^ ^ 

8s 



Us^ - hsM^ - 2M%s - M^) + sM^ In -^ 
6 M^. 



(C.6) 



and after a subsequent neglection of the proton mass 



TT 

48^ 



-s . 



(C.7) 



D Tables of integrals for the approach of section 3 

D.l First series of Zi(2)-iiitegrals 

This first series of integrals over zi or Z2 is used for the common treatment of the cross section 
in leptonic, mixed, and hadronic variables. The first integration is over zi or Z2- It corresponds 
to the angular integration over the angle if of the photon in the rest system of the proton as 
discussed in A.l. At this stage of the integration, there is no infrared problem. The necessary 
integrals are: 



[A. 
[1]. 



1 



^1(2) J 

1 

^1(2) 



TT 



,max 7 

-1(2) azi{2) 



nim 
1(2) 



^ 



A. 



^1,2 
1 

\/C'l(2) 

-^1(2) 

3/2 ■ 
(2) 



ct 



(D.l) 
(D.2) 

(D.3) 
(D.4) 



Using the relation Zi = Z2 + Q\ — Qf, one may derive 

1 

Z2 



Z1Z2 Qh — Qi 



1 

Zl 



(D.5) 



and thus reduce the corresponding integral type to known ones. The functions A, B, C are 
defined in (A.29)-(A.31). The integration boundaries zj^2) ™ may be found in (B.9). 

In appendix B.1.2, we perform the integration over zi(2) with dedicated cuts for the photonic 
energy and scattering angle. In this case, the above table of integrals becomes modified: 



[l]i 



^1(2) 
1 



Z 



1(2) 



1 



tta/ A, 



asm- 



5i 



(2) 



^1(2) Aq 



^d: 



I asm 



^l(2)-Bl(2) — C'i(2) 



ttVCi 

1 Y-Rz(^i(2) 

TtCi Zi(2) 



Zi{2)VDI 

l\ I ^1(2) 
'-1(2) 



^1(2) 



(D.6) 
(D.7) 
(D.8) 



95 



In the table, the abbreviation 



"max 
1(2) 

"min 
1(2) 



(D.9) 



is used, where the boundaries z"^/2)^'^^ have been derived in B.1.2. Further, we use (A. 28) in 
the shghtly rewritten form (B.8) and a similar relation holds for Dz'. 



D. 



AsW 



QT-'iW^ Ql Ml) - Q}\ \Ql - Qr\W\ Ql Ml) 



(D.IO) 



D.2 The integrals for the infrared problem 

For the angular integration, 

1 /•! 



[Ai 



2J~i 



d^A, 



(D.ll) 



one needs the following integrals which are calculated with account of the exact kinematics: 




— ln^ + ^' 



2P, 1-A' 
1 



21n2-2, 



ln2iln.l+''- 



+ 



A 1 - A 2A 



Li2 



' 2A . 



Li2 



' 2A . 

.A + 1' 



^f21n2-llni±^^ 



(D.12) 
(D.13) 

(D.14) 
(D.15) 
(D.16) 

(D.17) 



Since the definition of the velocities A depends on the lorentz system, some of the above 
integrals do so. 

Further, integrals over the Feynman parameter a are used: 



[A]. 



The first one is 



-A;2 



daA. 



ZL/jjj. 



with 



— A;^ = Q^a{l — a) + m^. 



(D.l^ 



(D.19) 



(D.20) 



and the Lm is defined in (4.29). Due to the invariance of —k"^, the integral (D.19) is relativistic 
invariant. 
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A second integral of the type (D.18) occurs in (4.23). It depends on the system in which 
it is calculated: 



I3a 



.U2 



-In 



(3o 



l+Pa 



g2 + 2m2 



iS$, 



^» - \i'^' -'-") i: ^- 



da 
I3a{-kl] 



In 



I -Pa 



(D.21) 
(D.22) 



The exact answer for this type of integrals may be found following the hints in [8]. For the 
proton rest system, the explicit result in the ultra-relativistic limit is (4.31). For the case of 
photoproduction, which was calculated in another rest system, the exact expression has been 
calculated numerically. Finally, we mention that for the case of the (76) rest frame, which 
was used e.g. in the case of Jaquet-Blondel variables, the kinematical simplifications are such 
that the S'$-function need not be introduced; see section 7.3. 

D.3 The second and third integrations for 5^^^ 

In this appendix, integrals over (first) Q^ and (then) y are collected. They are used for the 
calculation of the hard part of the infrared divergent cross section contribution. 

D.3.1 Leptonic variables 

The necessity to work out the absolute values in (B.39)-(B.40) explicitly for a subsequent 
integration over yh forces us to perform the integrations in the following four regions separately: 



gion I : 


„,min 
Vh 


< 


Vh 


< 


Vhl, 


Qhll 


< 


Ql 


< 


VhS, 


gion II : 


Vhl 


< 


Vh 


< 


yh2, 


Qhll 


< 


Q'i 


< 


VhS, 


gion III : 


yh2 


< 


Vh 


< 


Vhd, 


Qlii 


< 


Ql 


< 


VhS, 


gion IV : 


Vhd 


< 


Vh 


< 


C-(e), 


Qlii 


< 


Ql 


< 


Qli- 



We introduce now the operation 



{^}c 



S 






dQllAl 



where the operation [A]z is defined in (D.l). 
The integrals {vAjn^ i^i the regions I-IV are: 



(D.23) 



Ql 

Z1Z2 

IT 

ZIZ2 

IT 



Qz 



Qz 
II 



Qz 
II 



Qz 



1 ^^ b 

yih (-ai)(-a2)' 
2^1 



m 



0, 



Qz 




yih 



(D.24) 
(D.25) 
(D.26) 
(D.27) 
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m 

-4 



0, 



Qz 



Q 



Z\Zl 



m 



Q 



Qz 

J Qz 
\ IV 



^ b 



1 
.2 ^ III 



m 



Zn 



Qz 



1 



2:12:2 



li„« 



with 



I m 

ai 



2^ IV 



yih ^ 

.2^ IV 



Qz 



m 



zi 



Qz 



yih 



{l-yixi){yh-yhi), 
[l-yi{l -xi)]{yh-yh2), 
{yh-xivifil-yi -xiyir), 

yi - yh- 



Now we define the third integration over yh'- 



(A) 



yQz 



rVhi T fyh2 TT fVhd 

/ . dyn {AYq, + / dy^ [AYl^, + / dy^ {A} 



III 

Qz 



+ 



The result of this: 



^max(^) 



Vhd 



dyh {A}q^ 



.^^^^JyQz 



In^ln 



yf(i-x,Yii-m)s' 

m2 "^ 4M2e2(l - yixi)[l - yi{l - xi)] 

-yiil -xi) 



+ hn\i-m^ ^-^ 



2 
Li2 



In^ 



2 



yixi 



m 

A 



m 



yQz 



yQz 



In 



In 



1 -i//a;i)[l-yi(l -xi)] 
yi{l -Xi)S 



Li2(l), 



2Me(l - i/jxO ' 

yi{l -xi){l -yi)S 
2Me[l-yi{l-Xi)y 



From (D.35)-(D.37), (4.35) follows immediately. 



(D.28) 
(D.29) 
(D.30) 
(D.31) 
(D.32) 



(D.33) 



(D.34) 



(D.35) 
(D.36) 
(D.37) 



D.3.2 Leptonic variables. Photoproduction 

The first integration of (4.5) with Q^ = Qf, is over Z2- Then, one has to integrate (4.12), 
which becomes (5.4) for leptonic variables, over Ql within the integration region as shown in 
figure 36. We begin with these two integrations: 



{^}c 



i2 mill 



dQllAl, 



(D.38) 
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where the boundaries of the integration over Ql are (B.46), and the [A]z is defined and 
calculated in appendix D.l. Two such twofold integrals have to be performed: 



ior 



1 



Qz 



m 



-2T 



m 



W^-Mll3 
1 



-L/3, 



^^ ) Qz K '^2 ) Q2 

The L^ and (3 are defined in (6.8). There is only one third integration: 

f 1 1 _ r^h™ .,.o 1 

M 



W^-Ml 



^2min 



^'''iw^^M 



h 



In 



W'^-{M + m^) 



(D.39) 
(D.40) 



(D.41) 



The boundaries for this integration are (B.45). After three integrations, we finally arrive 
at (6.7). 



D.3.3 Mixed variables 

W consider the case Xm < 1- From the discussion of the roots of Ci and C2 in appendix B.3.2 
it may be seen that there are three integration regions: 



Region I : 


i/r(e) 


< yi 


< yid, 


Qlii 


< 


Ql 


< 


Qli. 


Region II : 


yid 


< yi 


< yi2, 


Qlu 


< 


Ql 


< 


yhS, 


Region III : 


yi2 


< yi 


< yT^, 


Qlii 


< 


Ql 


< 


yhS, 



with Ql,,, given by (B.27). 

Using the definition (D.23) for {Ajg^, one gets the following results in the three regions 



I-III: 



or 

Z1Z2 

IT, 

ZIZ2 



Qz 
I 



Qz 
II 



Qz 
II 



Z 



1 ) Qz 
III 



l-l 

[ Z1Z2 J 



m 



Qz 



z 



m 

-4 



1 ) Qz 
2^111 



yih "^^ 






Qz 



yih 



1 
yih 



2 In 



Q'l 



In 



[1 - x^{l + yh - yi)]{yir^ 



yi] 



m" 



. L Xr. 



yi - xiyir) 



m 

-4 



Qz 



1 

yih 

1 

yih 

0. 



..4 



1 

yih 

hln 



1 - Xm{l + yh -yi) 
yi - yi2 



(D.42) 
(D.43) 
(D.44) 
(D.45) 
(D.46) 
(D.47) 
(D.48) 
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The third integration is defined by 

J mm J,^ j j 



The results of the integrations are: 



^ZiZ2/ 



yQz 



m" 



4e2M2 



ln(l - yh) ln(l - Xr, 



m 



m 



yQz 



yQz 



In 
In 



{i-yh)S 

2eM ' 

(1 -yfc)(l -x„,)S 
2eM 



Xmjl -yh) 

Xm, -L 



From (D.50)-(D.52), (4.36) immediately follows. 



(D.49) 



(D.50) 
(D.51) 
(D.52) 



D.3.4 Hadronic variables 

For the hadronic variables one has to integrate 5^^^ over the kinematical region of figure 41. 
With the aid of the boundary conditions which are derived in appendix B.4.2 the integration 
region may be split into the following three subregions: 



gion I : 


i/r°(e) 


< 


yi 


< 


yi3, 


Qh 


< 


Ql 


< Q?i, 


gion II : 


l/«3 


< 


yi 


< 


yi2, 


Qiw 


< 


Qi 


< QU. 


gion III : 


yi2 


< 


yi 


< 


Vh 


Qiw 


< 


Qi 


< Qini 



The operation {-AJq^ is defined as follows: 

{A}qz = S J-^^dQMl 

The integrals {A} in the three regions I-III are: 



Q 



I 



^1^2 J Q 
( 2^1 

m 



z 



m 

~4 



yih ^ 
1 



yihi'^-yih) 

i-yi 
yihi^-yh)' 



\ziZ2 



II 



1 

yih 



21n^ + ln 



{l-yi){2 + ^,/S-yh) 



m^ 



[2(1 - yi) + ^/Xh/S + yh]{l -yh- Xhyhr) 



m 



yihi^-yihY 



(D.53) 



(D.54) 
(D.55) 
(D.56) 
(D.57) 
(D.58) 
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771 

'4 



II 



i-yi 
yihC^ -VhY 



I ^1^2 Jo. 



771 



III 



1 

yih 

0, 



In^ + ln 



{^-yif 



771^ 



XhVhr- {I - yi)il - yih) 



771 



2^111 



1 -yi 



[4}q, yihi^-Vh)' 

The third integration is: 



f'Vl rVl f'Vl 

(^),Q. = 1, , dy,{AYQ, + / ' dyAA}l, + / dy,{AYS.. 



And the resuhs of the integrations are 



Q 



'A 



i.«i„(^-«"'^"^ 



i2\ 



m 



771 



In 



j/Qz 



j/Qz 



m- 4e2M2 

2eM y^ r 



2eM 



Li2(l-t/0-Li2(l), 



From (D.64)-(D.66), (4.37) follows immediately. 



(D.59) 
(D.60) 
(D.61) 
(D.62) 



(D.63) 



(D.64) 
(D.65) 
(D.66) 



D.4 The second integration for 5^. Mixed and hadronic variables 

When performing the integration over yi at given values of y^-, Q\i and Qj within the limits 
ymm g^j^^ yi"^"-^, one faces eleven types of integrals which should be treated separately in 
regions I and II of figure 39. 

For details of their definitions see the appendices B.3.3 for the case of mixed and B.4.3 for 
hadronic variables. 



The integrals are: 



[Ai 



S 



y'l" 



Vl 



Integrals for region I with Qf > Qf^ 

I 



2. 



In 



Q 



^VhQ 



2 ' 



^ 



4^^ 



Qf 



1- 






QJ 



771^ 



Ql 
Of 



dyiA. 



(D.67) 



(D.68) 



(D.69) 
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^ 



yi 

/cl 

yi 



m 



2 



-ii 

- y 
I 



y 



-'y 
I 

y 
I 

y 
I 



4'" 



1 -yfc 

Of 



Q 



yhQl 



Qf 
Ql - yhQ^ 






+ 



Q1-QI + yhQl 



Q} 



I ^Q 



Ql + yhQl 



.3/2 



Cl 



m 

3/2 



Qi 

S\Ql-ynQl? 

S\Ql-yhQl)\ 
Ql 



Ql 



/c; 



nl 



-ij/ 



+ 



Ql 



Ql. 



'Ql 

.^Qt 
yhQl Ql 1^ 
2Q^ Q^ 2, 



252(g2 

0, 



i//.Qf)(Q/ 



Ql)^' 



y 



cr 

LCf 



Q1-QI + ynQl 



-'y 



2S\Ql 



0. 



yhQDiQi 



Qir 



Integrals for region II with Ql < Q\: 
-. nil 



In 



g-i y 



II 



.yh, 



/CT 



yi 

/C'o 



II 

-I J/ 
II 




Q' 



yh ^Qi 



qI 

' Ql + yhQl 



Ql 
:i-yh){2Q 



Ql 



/CT 



Ql) , Ql- Ql + yhQl 



-'y 
■ill 

y 
■ill 



Qt 
:i-yk)s' 

;i - y^rs' 



Ql 



Ql- 

' Ql 
Ql 



'^ + yK)Ql 



Q1 + 



Ql 



^1 



(D.70) 

(D.71) 
(D.72) 
(D.73) 
(D.74) 

(D.75) 
(D.76) 
(D.77) 
(D.78) 



(D.79) 



(D.80) 



(D.81) 







(D.82) 


1 

VC2J 


II 

y 


(D.83) 




(D.84) 






(D.85) 



102 



m 



II 



ct 



3/2 



0, 



m 



cT 



vfi^Vi 



CI 



■3/2 








1 


y 


2{l -y^){Qj-QlYS^' 


II 




= 


-- 0, 


y 




ii 


QI-Q1 + VhQl 



-'y 



2Ql{l-y,){Qf-QiyS^- 



(D.86) 
(D.87) 
(D.88) 
(D.89) 



E Tables of integrals for section 7 

E.l The cos i?i?, (/9i?- integrals 
The integrals are defined as follows: 



[A 


fR 


- 2J0 '^ 


^I'Pr. ~ 


= 1, 


[1-^4, = 


= A, 


\i-mr\ 


•fiR 


- ^'^>'- 



The following abbreviations are used: 
1 



A 
B 



Z2 



(l-yh){r + m') + 
It I v^i 



Z2 



Ql - Z2 - VhiQl + T + m^) 



COS-i?/ 



2tS 



A^sin-i^pSin-i^jj. 



(E.l) 

(E.2) 
(E.3) 

(E.4) 



(E.5) 
(E.6) 



After applying these exact formulae, in a second step the following integrals are applied: 

2 7-1 

1, 






A, 



[cos-i? 



^J1?iJ 



1 

■ Zl 

1 

~J2 



0, 
r 1 



(E.7) 

(E.8) 
(E.9) 



In 



Ql^T^rn^ ^ ^f\x 



Z2 \f\ '" Ql + T + m^ 
r 1 



Z2 m^Z2 



(E.IO) 
(E.ll) 



In ultra-relativistic approximation: 



, Qi + T + m' + ^i 

In -:z7, 7== ~ Lit-i 



Ql + r + m?- VAT 



(E.12) 
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and 



■1- 

-Zl. 



Ql 



L,. 



(E.13) 



The L^ and Q^ are defined in (7.16) and (7.17), respectively. 

The twofold integrals over the two photon angles, J7'^[v4.] = (r — 'm^)/(47rr) ^ (KIrA., with 
A = 1/zi and A = 1/zf agree exactly with those which have been calculated for bremsstrah- 
lung problems in e~^e~ annihilation in [80]. 

E.2 The third integration for 6^^.^ 

In the integration of the infrared divergent part of the cross section we use the following 
regularized integrals for the integration of the hard, finite contributions in hadronic variables: 



[A]'. 



dZ2A, 



(E.14) 



Z2 



with Z2 = 2me, and r"'"''' = Ql{l - yh)/yh, see (B.94). 

In the case of Jaquet-Blondel variables, the upper boundary differs: 



[A] 



JB 



dzoA, 



(E.15) 



with r"'^^ = (1 - a;jB)(l - ym)S, see (B.106). 
The integrals are: 






In 



2me yh , 



2i„«|,J«Li^Ui„.foii^' 



Z2 



■ ]^ 1 JB 

-Z2\z 



m" 



2me yh 



2Lio 



Li.rii 



m^ yh 



Vh, 



In 



1 -a;jB)(l-yjB)5' 



Z2 



2 In 



Q 



JB 



m^ 



In 



2mt 
^l-xjB)(l-yjB)5' 1 2 (1 -a;jB)(l -yjB)5' 



- 2Li, (1 - X. 



JBJ 



2rae 



-In^ 
2 



m^ 



where 



X 



JB 



1-xjb(1 -^jb) 



The regular integrals may be found in appendix E.3. 



(E.16) 

1), 
(E.17) 

(E.18) 
(E.19) 
(E.20) 



E.3 The third integration for 5^1 

A second series of integrals over Z2 is used in order to perform analytically even the last, third 
integration in case of hadronic variables. Here, this integration corresponds to an integral over 
the invariant mass r; see the discussion in section 7.4. 
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The finite integrals over Z2 or, equivalently, over r, are defined as follows: 

Vt = I drA, 



[At 



JB 



drA. 



(E.21) 

(E.22) 



Here, at the upper integration limit the ultra-relativist approximation may be applied: 



^max _ q2 ^x _ yi^)lyy^ in the case of hadronic variables, and r'' 
in the case of Jaquet-Blondel variables. 
The integrals are: 



[1]; 



Q 



2 ^-yh 

Vh 



-,h 



In 



-,h 



T 
1 

1 

Of 
1 

1 
Of 



[L.]; 



m"^ Vh , 



J T 

h 



Qi 






Vh 

^|21n- + (l-y.) 
Vh [ Vh 

In^y/j-lny/jln— I + Lig 



21n 



Ql 



In 



'Q^i-y/.^ 



m^ 



m" yh 



.Vh. 



-UM] 



J T 



Q^ 



2Qt 



-(l-2t/;,)lni/^ + (l-y;,) (inl 
-{l-2yl)\nyn + {l-yl)\4 



^m^l-yh 

Ql VH \ 
m? 1 -yh) 



{l-XiB){l-yiB)S 

(E.23) 
(E.24) 
(E.25) 
(E.26) 

(E.27) 
(E.28) 
(E.29) 
(E.30) 
(E.31) 



)+2 

+ {l-yK){2 + yh) 



Further, 





rl" 


JB 








= 


= In 




W _ 


r 




ri' 1 


JB 




I 


Tr. 


r 


= (1 



1 -a;jB)(l -yi^)S 



m^ 



ym) 



In'XjB + lnXjBln 



Q 



2 

JB 



"^^(1 -I/jb) 



Li2(X, 



JBj 



Li2(l) 



where the Xjb is defined in (E.20). 
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F Leading logarithmic approximations 

In the leading logarithmic approximation, the 0{a) corrections consist of initial state radia- 
tion, final state radiation, and the Compton peak contribution. 
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where the latter contribution exists only for leptonic and mixed variables. For leptonic vari- 
ables it is numerically important: 
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The other two cross sections have the following structure: 
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For the different sets of variables, the definitions of the x, y, S, as well as the lower integration 
boundary zq, differ. For leptonic variables, they are well-known since long; see, e.g. [17, 77]. 
For the cases of mixed and Jaquet-Blondel variables, they have been derived in [74] and for 
further sets of variables in [76] where also the formulae for higher order LLA corrections with 
soft photon exponentiation are collected and numerically discussed. 

The formulae which were used for our comparisons in section 8.2 may be found in table 1. 

In the case of Jaquet-Blondel variables, the hadronic final state is treated totally inclusive. 
Thus, in accordance with the Kinoshita-Lee-Nauenberg theorem [81] there is no LLA correc- 
tion from final state radiation. One should also mention that a gauge invariant separation of 
initial and final state radiation is possible for the leading logarithmic corrections but not for 
the complete order 0{a). 

In (F.3), it is 



X 



J{x,y,Q^) 

x{zo) 



yS' 
d{x,y) 



d{x,y) 



1. 



(F.4) 

(F.5) 
(F.6) 



From the last equation, the lower integration boundary derives. 

More details on LLA corrections in a large variety of variables may be found in [78] . 
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Initial state radiation 


Final state radiation 


Variables 


leptonic 


Jaquet-Blondel 


mixed 


leptonic 


Jaquet-Blondel 


mixed 


y 


z + y-l 

2, 


y 

z 


y 

z 


z + y-l 

z 


- 


y 


S 


zS 


zS 


zS 


S 


- 


S 


Q' 


zQ^ 


T^<?' 


zQ^ 


z 


- 


z 



Table 1: The definition of scaling variables for the different leading logarithmic cross sections. 
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